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HISTORY ANALYSIS OF MID-RISE STEEL-STRUCTURE BUILDINGS 
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Saeed Amiri 7®, Aram Saaed'®, Ali Yahyapour'® 
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2 Polytechnique Montreal, Department of Civil, Geological and Mining Engineering, Montreal, Canada 


Abstract: Response history analysis using a time integration method isa Aram Soroushian* 
powerful versatile tool in accessing structures seismic behaviours. To — B-mail: a.soroushian @iiees.ac.ir 
reduce the analysis run-time, a technique was proposed in 2008 for time 

integration with steps larger than the steps of ground motions. The 

technique has been implemented in seismic assessment of frames, Received: 11.01.2023 
buildings, bridges, silos, etc., leading to considerable reductions in the Revised: 30.01.2023 
analysis run-time, without notable effect on the response accuracy. The Accepted: 20.02.2023 
technique has recently been named as the SEB THAAT (Step- 

Enlargement-Based Time-History-Analysis-Acceleration-Technique). To 

use the SEB THAAT, the smallest dominant period of the response needs © The Author(s) 2023 

to be available prior to the analysis. In this paper, concentrating on 5-20- 
floor steel-structure buildings, a simple engineering comment is proposed 
that eliminates this need. As a result, in response history analysis of mid- 


rise steel-structure buildings subjected to ground motion, by using the This work is licensed under a Creative 
proposed comment, we may reduce the analysis run-time, significantly, Commons Attribution-NonCommercial 
without any initial information about the response. The reduction is 50% 4.0 International License 


for the linear analyses. 

Keywords: response history analysis, ground motion, mid-rise steel- 
structure buildings, the SEB THAAT, integration step, excitation step, 
run-time, accuracy. 


Introduction 


Response history analysis using a time integration method is a powerful versatile tool, for analysing 
structures subjected to ground motions, irreplaceable in important analyses [1-3]. After discretization in space, 
the initial value problem, representing the behaviour of the structural system, is expressible as [4-7]: 


Mii (t) + fine(t) = (MP iig(¢)) OS' Sta 


u(t = 0) = Ug 
Initial Conditions: |U(t = 0) = Ug : (1) 
fint(t = 0) = finto 


Additional Constraints: Q 


In Eq. (1), t and teng imply the time and the analysis time interval; M is the mass matrix; fj,; is the vector of 
the internal force (in linear problems, generally fj,;_ = Ku + Cu, where K and C stand for the matrices of 
stiffness and viscous damping, respectively); ti,(t) implies the single-component ground acceleration, and 
is a vector with the size of the degrees of freedom, needed for matrix multiplication and considering spatial 
changes of i,(t); u(t); u(t), and ii(t), denote the vectors of displacement, velocity, and acceleration, 
relative to the ground; Uo, Uo, and fint,, define the initial status (generally all zero), and Q implies the 
limiting conditions due to nonlinearity. (When the ground motion is multi-component, I’ and iig(t) 


will change to a matrix and a vector, respectively). The core of response history analysis is an approximate 
step-by-step computation, widely known as direct time integration (see Fig. 1 and [1,3,6]). The resulting 
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Fig. 1. A pictorial review of direct time integration analysis of Eq. (1) [3] 


responses are inexact, and the analysis run-times are generally considerable [3,6]. The integration step is the 
analysis parameter, with adverse effects on the run-time and the accuracy [3]. Accordingly, a regulation for 
assigning an appropriate value to At is essential. A simple comment, broadly accepted in practice, is as follows 
[2,3,8]: 


T 
At =< Mi (5. h, at) : 2, 
inly f (2) 


In Eq. (2), T is the smallest dominant period in the time history of the response, h implies the largest integration 
step preserving numerical stability and consistency, ,At denotes the step, by which, the ground motion is 
digitized, and X equals 10 for linear systems, 100 for nonlinear systems not involved in impact, and 1000 for 
systems involved in impact. When, the dominant term in the right hand side of Eq. (2) is At, the difference 
between ,At and the next smallest term in the right hand side of Eq. (2) implies the effort required mainly to 
account for all the excitation data. In order to eliminate or lessen this effort, in 2008, a technique was proposed 
[3,9,10], that replaces the earthquake record with a record digitized in larger steps. The technique, which is 
recently named “the SEB THAAT (Step-Enlargement-Based Time-History-Analysis-Acceleration- 
Technique)” [10], is formulated such that to prevent any negative effect because of the record replacement on 
the convergence of the computed response to the exact response. Considering that convergence is the main 
essentiality of approximate computations [11,12], the expectation from the SEB THAAT is to reduce the 
computational effort with negligible effect on the response accuracy. In view of the studies carried out since 
2008, the SEB THAAT has been successful, in earthquake engineering and beyond; see Table 1 [3,10] and 
[13]. Nevertheless, to set the scaling value of the record’s step (i.e. how much to enlarge the step), some 
information about the response is needed, prior to the analysis [3,10]. The objective in this paper is to eliminate 
this need for a special class of analyses. In view of the number and social importance of buildings with 5-20 
floors, hereafter referred to as mid-rise buildings, seismic response history analysis of steel-structure mid-rise 
buildings is considered as the special class of the analyses. Accordingly, by achieving the objective, the 
computational effort needed for a large number of response history analyses will be simply reduced. The 
simplicity in increasing the efficiency and the amount of the increase can encourage engineers to use response 
history analysis in practice, as well. For achieving the objective, the main attention is paid to Eq. (2), the current 
design and analysis practice, and the trends of the future advancements. A brief review on the SEB THAAT 
is presented first. An engineering comment for selecting the enlargement scale is introduced next. 
The effectiveness of the comment is tested afterwards, and finally, the paper is concluded with a discussion on 
practical issues, as well as an overview of the achievements. 
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Table 1. Some past tests on the SEB THAAT with regard to seismic response history analysis 


System Reduction in run-time (%) 
A thirty-storey building 50 
3-component earthquakes 66 

Silo 77 

Water tanks and Silos 66 

Bridges 45-80 
Residential buildings 50-87 

An earth dam < 80 

Milad telecommunication tower 50-70 

Different lifelines 50-90 


A brief look at the SEB THAAT 


An overview of the SEB THAAT, its formulation, limitations, challenges, and future prospects, is recently 
presented in [10]. Therefore, for brevity, only application of the SEB THAAT to direct time integration 
analysis is referred to in Fig. 2, and the most important features are reviewed as follows: 


1. Define the structural model and the excitation (the excitation digitization step equals -At) 

2. Select the integration method 

3. Select the details of the nonlinear solution (if needed) — — — Direct time integration 

4. Select the integration step At, regardless of the SEB THAAT Direct time integration after 
5. Step-by-step direct time integration using At as the integration step applying the SEB THAAT 
6. Assign a value to the enlargement scale n 

7. Use the SEB THAAT to change the excitation to an excitation digitized at step n At 

8. Step-by-step direct time integration using n ,At as the integration step 


Fig. 2. Application of the SEB THAAT to arbitrary direct time integration analysis 


1. The SEB THAAT has a mathematical basis, with the aim and formulation to preserve the responses’ 
convergence, without necessarily avoiding changes in characteristics of the earthquake record (see [9]). 

2. Considering Eq. (2) as a basis for selection of the integration step, the SEB THAAT is to be applied, only 
when ,sAt dominates the right hand side of Eq. (2), i.e. 

At < Min (. h) (3) 

It should however be noted that Eq. (2) is not rigorous [2,3,8,14], and hence, the SEB THAAT may be 
successful, even when Eq. (3) is invalid; see for instance [15]. 

3. The formulation of the SEB THAAT enables consideration of arbitrary real number greater than one as 
the step enlargement scale [16]. 

4. With appropriate details of nonlinear solution and sufficient computational facility, the SEB THAAT can 
be as successful in application to nonlinear analyses, as it is in application to linear analyses; see [17]. 

5. While in linear analyses application of the SEB THAAT causes reduction in the analysis run-time, R, 


obtainable from 


n-1 
R = 100 


% (4) 


(n is the step enlargement scale), in nonlinear analyses, the reductions are not necessarily obtainable 
from Eq. (4) [3]. 
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6. With a basis on convergence of the computed results, the SEB THAAT may be successfully applicable in 
fields different from structural dynamics and earthquake engineering [13]. 

7. Compared to the SEB THAAT, direct down sampling [18] is less effective, i.e. the accuracy of the target 
response obtained from the analysis when using the SEB THAAT is more than when using direct down 
sampling instead of the SEB THAAT [19]. 

It is also worth noting that, the future of the SEB THAAT is promising, due to its simplicity, remarkable 
effectiveness, and everyday more availability of different data as digitized records with smaller digitization 
steps. 


A new engineering comment 


For the SEB THAAT to reduce the computational effort, a positive value larger than one should be assigned 
to the step enlargement scale, n; see also [3,10]. Besides, the greater the n, the more will be the reduction in 
the run-time, especially for linear analyses [3,9,10]; see Eq. (4). Accordingly and in view of Eq. (2), for the 
analysis most efficiency and good accuracy of the response, the excitation step n ¢At_ should govern the right 
hand side of Eq. (2). Accordingly, when h > © (broadly recommended [3,4,6,20,21]): 

natal 
n= X;At’ (5) 
which is effective, when the resulting n is greater than one. (The reason of using an approximation sign in 
Eq. (5) is that, different from h and At , the T/X in Eq. (2), the definition of T, the values of X, and the 
form of Eq. (2), are not rigorous (see [3,14])). 

The existing seismological instrumentations provide the capability of recording ground motions, in steps, 
as small as 0.004 sec [22]. These instrumentations are in continuous progress, towards smaller digitization 
steps; see [23]. Considering this, along with the about largest digitization steps currently in use [24], we can 
conclude that: 


At < 0.02 sec, (6) 
and, 
ale Pa 207, . (7) 
X At X 


Therefore, considering that smaller values of n will be more reliable from the standpoint of response accuracy, 

it is practically acceptable to replace Eq. (5), with 

50T 

—. (8) 
X 

Accordingly, if for a class of structural analyses there exists a minimum for T (see Eq. (2)), i.e. Tynin » Such 

that: 


n= 


20 Tmin 
= 
the n is suitable in application of the SEB THAAT in that class without any information about the response. 
Taking into account 5-20-floor steel-structure buildings, designed according to the Iranian codes!”, and in 
agreement with the International code? and [25,26], the least dominant periods of the displacements linear 
oscillations, because of ground motions, satisfy (displacements, velocities, and accelerations, are the unknowns 


n Si. (9) 


generally being computed first in time integration analysis): 


Tmin = > 0.38 —0.40sec . (10) 


' BHRC (Building and Housing Research Centre), Standard No. 2800-05, Iranian Code of Practice for Seismic Resistant 
Design, Iran, 2007 (in Persian). 

2 INBR (Iranian National Building Regulations), Iranian National Building Code, Part 10-Steel Structures, Iran, 1993 
(in Persian). 

3 ICC (International Code Council), IBC - International Building Code, Club Hills, USA, 2003. 
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Besides, nowadays the technology of buildings construction proceeds towards lighter designs, and generally 
less stiff structural systems and larger values of Tyjin [27]. Considering these, in application of the SEB 
THAAT to analysis of 5-20-floor steel-structure buildings, the following comment is reasonable: 

w= 2: (1) 
It is worth noting that Eq. (11) is in agreement with the experiences on the SEB THAAT’s application, 
reviewed in [10]; see also Table 1. Equation (11) implies an engineering comment, for implementation of the 
SEB THAAT [9] which, is considered above for linear analyses, and simply reduces the analysis run-time for 
about 50%. Nonlinearity because of inelastic behaviour is inherent to structures subjected to severe earthquakes 
[1], and generally causes increase of both T,,;, and X. Besides, the T/X in Eq. (2), the definition of T, the 
form of Eq. (2), and especially the values of X are not rigorous [14]. Furthermore, seismic standards 
mostly permit linear analyses for ordinary building structures***7*-'°, and the only seismic standard, with a 
procedure for nonlinear response history analysis, i.e. NZS 1170.5:2004"', also proposes a regulation to check 
the accuracy of the response. Accordingly, it is reasonable to use Eq. (11) in application of the SEB THAAT 
to analysis of 5-20-floor steel-structure buildings inelastic behaviours, as well. Nevertheless, because of the 
iterative nonlinear solutions [7,28,29], using Eq. (11) will not lead to twice faster analysis when the behaviour 
is inelastic. Consequently, we can use Eq. (11) in application of the SEB THAAT to seismic response history 
analysis of 5-20-floor steel-structure buildings, and reduce the analysis run-time, without prior information 
about the response. This claim is tested in the next section, considering quiescent condition at t = 0 and S.L 
system of units for all computations. 


Numerical study 

Simple test on a six-floor steel shear frame 

Consider response history analysis of the system defined in Fig. 3 (g stands for the acceleration of gravity) 
and Table 2, by the C-H method [30] (p,. = 0.8). The top displacement is set as the target response. The 
suitability of using Eq. (11) is tested via analysis with the integration step At = ,At and a second analysis 
using At = 2 ,At . The sufficient accuracy of the response and the validity of Eqs. (3) and (10) are displayed 


in Figs. 3 and 4, while the numbers of integration steps (37500 and 18750, for Figs. 4(a) and 4(b), respectively) 
imply 50% reduction in the analysis run-time. 


Floor 6 


Floor 5 Shear Frame 


oa 


Floor 4 


pAt=0.004 sec 


Floor 3 


Floor 2 


Floor 1 


0 30 60 90 120 150 


(a) <*> (b) Time (sec) 


Fig. 3. Pictorial introduction to the first example: (a) Structural model, (b) Ground motion 


4 NZS (New Zealand Standards), NZS 1170.5: 2004 Structural Design Actions, Part 5: Earthquake actions, New Zealand, 2004. 

> BHRC (Building and Housing Research Centre), Standard No. 2800-05, Iranian Code of Practice for Seismic Resistant 
Design, Iran, 2007 (in Persian). 

® INBR (Iranian National Building Regulations), Iranian National Building Code, Part 10-Steel Structures, Iran, 1993 
(in Persian). 

7ICC (international Code Council), IBC - International Building Code, Club Hills, USA, 2003. 

8 BCJ (Building Centre of Japan), Structural Provisions for Building Structures, Tokyo, Japan, 2001. 

° EAK, Greek Code for Seismic Resistant Structures, Athens, Greece, 2000. 

‘0 NRCC (National Research Council Canada), National Building Code of Canada, Canada, 2005. 

'! NZS (New Zealand Standards), NZS 1170.5: 2004 Structural Design Actions, Part 5: Earthquake actions, New Zealand, 2004. 
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Table 2. Main properties of the structural model in the first example 
Floor (i) / Mode (i) 
Property 
1 2 3 4 5 6 
Mass ( mm, )x10” 1.8 1.8 1.8 1.8 0.6 0.6 
Stiffness (k,)x 10! 1.20 1.20 1.20 1.20 0.20 0.20 
Natural period (7. ) 2.7191 1.4919 0.7695 0.6544 0.4973 0.4088 
Damping Classical damping [2], considering 2% damping for the 1“ and 3" natural modes 
0.6 0. 

q 0.2 oe | | my q 0.2 \ iv 13 

F | Th, 5 lhw } 

Q se tihiL a att * { g ad pg 4H : * " { 

=, nie | 2 Ty 1 || {I 

g 0.2 mahi us 4 02 Hil "Wy! 

a ‘he a | | 

a0 |} a 04 1} 

= |! ea | | 

0. 4 0.6 4 
0 x 60 90 120 150 x0 60 90 120 150 
(a) Time (sec) (b) Tume (sec) 


Fig. 4. Target response computed for the first example by the C-H time integration method (p.. = 0.8) 
using an integration step equal to: (a) fAt , (b) 2 At (by means of the SEB THAAT) 


A fifteen-floor steel-structure building 


Consider the 15-floor steel-structure building displayed in Fig. 5. There are two identical axes of symmetry 
in the plan; the lengths of the spans are four meters; the bracings are placed at the last spans of the surrounding 
frames; the structural system is dual!?, and the other basic details are reviewed in Table 3. The structure is 
designed in a previous study [31], based on the Iranian standards'*""*, for the lowest total cost of construction 
using a fully constrained optimal criterion [31,32]. The members’ cross-sections are as reported in Table 4. 

Four two-component records, of the historically most devastating earthquakes in Iran, are considered, as 
the ground motion records (see Fig. 6) (for the first earthquake, the two components are identical). 
The structure is modelled as a three dimensional shear frame [2,24]. The translational natural periods are 
in the interval (0.035 sec 0.8 sec). The seismic response history analyses are carried out using the average 
acceleration time integration method [33], twice, for each two-component record; once, ordinarily, and once, 
after application of the SEB THAAT [9,10], considering Eq. (11). The acceleration at top, the mid-height 
displacement, and the base shear, are considered as the target responses. The computed histories are depicted 
in Figs. 7-10, where, for further clarity, the mid-sections of the time histories are not displayed in Figs. 8-10. 
These figures show that Eq. (11) may present an adequate comment for assigning a value to n in application 
of the SEB THAAT in response history analysis of mid-rise steel-structure buildings. 


2 BHRC (Building and Housing Research Centre), Standard No. 2800-05, Iranian Code of Practice for Seismic Resistant 
Design, Iran, 2007 (in Persian). 

3 Tbid. 

'4 INBR (Iranian National Building Regulations), Iranian National Building Code, Part 10-Steel Structures, Iran, 1993 
(in Persian). 
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(a) (b) 
Fig. 5. The building structure in the second example: (a) Schematic view, (b) Top view 
Table 3. Basic details for the structural systems in the second and third examples 
Material Steel (ST-37) 
Occupancy Residential (in Shiraz, Iran; mass: 175000 and 200000 Kg for the roof and other floors respectively) 
Seismic zone factor | 0.35! 
Soil type II (375 m/s < Vs <750 m/s )'°# 
Floor height 3 meters 
Bracing x 
Damping Negligible (considered zero in the analyses) 


“ Vz is the velocity of shear waves. 


Table 4. Members’ cross-sections for the structural system introduced in Fig. 5 and Table 3 


Floor Inner Peripheral Corner Inner beams Peripheral beams Bracings 
(from columns columns columns 
ground) 

1 IPB400 * Box400*12.5 | Box360*10 | IPE 300 2IPE160 2L130* 12 
2-3 | IPB400? Box400*12.5 | Box360*10 | PL200*10 + 2PL240*20 2TPE200 2L130*12 
4-5 |Box320*10 | Box400*16 | Box360*10 | 2IPB200 2IPE200 L180*16 
6-7 | Box320*10 | Box360*16 | Box320*10 | PL200*10 + 2PL240*20 2TPE200 2L110*10 
8-11 | Box320*10 | Box360*12.5 | Box320*10 | PL200*10 + 2PL240*20 INP240 2L110*10 

12-15 | Box320*12.5 | Box320*10 | Box280*8 | PL150*10 + 2PL180*15 PL220*6 + 2PL120*10 | L130*12 


* For construction problems, the properties in the weaker direction of the cross-section are considered as the properties in both directions. 


'S BHRC (Building and Housing Research Centre), Standard No. 2800-05, Iranian Code of Practice for Seismic Resistant 
Design, Iran, 2007 (in Persian). 
'6 Tbid. 
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Fig. 6. Two-component records of the four most devastating earthquakes in Iran, as the ground motion 
records in the second example: (a) Naghan (1977), (b) Tabas (1978), (c) Abbar (1990), (d) Bam (2003) 


Sixty-five 10-20-floor steel-structure buildings 


In this section, the adequacy of Eq. (11) is studied, in view of 25 ten-, 25 fifteen-, and 15 twenty-story steel- 
structure buildings, each with two identical axes of symmetry, X bracings with the configuration displayed in 
Fig. 11, and span lengths constant throughout the structure, equal to either of the followings (see also [31]): 

s=4, 5, 6 (meters). (12) 
Two two-component ground motions, selected, based on the soil type and shear wave speed, are applied, at 
the ground level, in the principle directions of the structures (see also Fig. 12 and Tables 3, 5, 6). The response 
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Fig. 8. Target responses computed for the second example when subjected to the records in Fig. 6(b): 
(a) The starting ten seconds, (b) The ending ten seconds 


history analyses are carried out using the average acceleration method [33], once with the step At = ,At, and 
then with the step At = 2 At (after applying the SEB THAAT considering Eq. (11)). The maximum relative 
difference between the two responses in the L,,. norm [34] is reported in Tables 7-9. Smallness of the reported 
values, that the ordinarily computed responses are not exact [2-4,6,35], and the 50% reduction in the analysis 
run-time, imply the good performance of the SEB THAAT when using Eq. (11). 
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Fig. 9. Target responses computed for the second example when subjected to the records in Fig. 6(c): 
(a) The starting fifteen seconds, (b) The ending fifteen seconds 
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Fig. 10. Target responses computed for the second example when subjected to the records in Fig. 6(d): 


(a) The starting ten seconds, (b) The ending ten seconds 
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Fig. 12. Records of the ground motions in the third example: (a) Manjil (1990), (b) Koujoor (2004) 
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Table 5. Systems of the buildings structures in the third example according to the Iranian codes'”'8 * 
System 10-floors | 15-floors | 20-floors 
Dual consisted of special moment frames and concentrically braced frames + + + 
Special moment frames + + + 
Dual consisted of intermediate moment frames and concentrically braced frames + + + 
Intermediate moment frames + + - 
Ordinary concentrically braced frames + + - 


a “+” and “-” imply being and not being considered in the design code (and in this paper), respectively 


Table 6. Groups of identical structural members in the sixty-five buildings in the third example 


Buildings 


Floors with identical structural members *° 


10-floor buildings 1-2-(3, 4)-(5, 6)-(7, 8, 9, 10) 


15-floor buildings 1-(2, 3)-(4, 5)-(6, 7)-(8, 9, 10, 11)-(12, 13, 14, 15) 


20-floor buildings 


1-2-(3, 4)-(5, 6)-(7, 8)-Q, 10, 11, 12)-(13, 14, 15, 16)-(17, 18, 19, 20) 


a The numbers in each “(_)” address the floors’ numbers with identical structural members (groups in height) 
b The seven groups with identical structural members in plan are: internal columns, side columns, comer columns, internal 
beams, peripheral beams, internal bracings (only for cases in Figs. 11(b) and 11(c)) and side bracings 


Table 7. Differences in the L., norm, between the responses obtained when using the SEB THAAT and 


Eq. (11) and the responses computed using At = ,At, for the 10-floor buildings in the third example (%) 


Structural system Top Mid-height Base 

acceleration | displacement shear 
Dual consisted of special moment frames and concentrically braced frames £5.39 < 0.83 < 1.79 
Special moment frames < 5.92 < 0.89 < 0.72 
Dual consisted of intermediate moment frames and concentrically braced frames < 5.46 < 0.66 < 1.91 
Intermediate moment frames < 5.77 S 0.83 $1.17 
Ordinary concentrically braced frames < 6.20 < 0.89 < 2.55 


Table 8. Differences in the L,, norm, between the responses obtained when using the SEB THAAT and 


Eq. (11) and the responses computed using At = At, for the 15-floor buildings in the third example (%) 


Structural system Top Mid-height Base 

acceleration | displacement shear 
Dual consisted of special moment frames and concentrically braced frames £0.53 < 1.67 < 5.98 
Special moment frames < 0.15 < 0.57 < 5.43 
Dual consisted of intermediate moment frames and concentrically braced frames < 0.60 < 1.80 < 5.96 
Intermediate moment frames < 0.13 < 0.59 < 5.90 
Ordinary concentrically braced frames < 0.62 < 1.46 < 5.58 


Table 9. Differences in the L,, norm, between the responses obtained when using the SEB THAAT and 


Eq. (11) and the responses computed using At = ,At, for the 20-floor buildings in the third example (%) 


Structural system Top Mid-height Base 

acceleration | displacement shear 
Dual consisted of special moment frames and concentrically braced frames £0.23 < 0.79 < 5.38 
Special moment frames < 0.14 < 0.98 < 5.63 
Dual consisted of intermediate moment frames and concentrically braced frames < 0.28 < 1.15 < 5.13 


'7 BHRC (Building and Housing Research Centre), Standard No. 2800-05, Iranian Code of Practice for Seismic Resistant 


Design, Iran, 2007 (in Persian). 


'8 INBR (Iranian National Building Regulations), Iranian National Building Code, Part 10-Steel Structures, Iran, 1993 


(in Persian). 
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A 15-floor steel-structure building with nonlinear behaviour and non-classical damping 


A 15-floor steel-structure building model is under study for the top displacement and base shear; see 
Fig. 13 and Table 10. The behaviour is nonlinear, with attention to Fig. 14, with regard to which, for the top 
displacement, T = 1.0 > 0.40 sec; see also Eq. (10). This is also in agreement with the fact that the first seven 
natural periods of the system are greater than 0.4 sec (see the mode shape in Fig. 15), nonlinear behaviour can 
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Fig. 13. Pictorial introduction to the fourth example: (a) Structural model, (b) Ground motion 
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Fig. 14. Exact responses of the fourth example 
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Fig. 15. Shape of the seventh natural mode (Tz = 0.42 sec), for the linear system 
corresponding to the nonlinear system introduced in Fig. 13(a) and Table 10 
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Table 10. Main properties of the structural model in the fourth example 
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be considered an extension of linear behaviour, the first few natural modes play the main role in the structural 
behaviour, and that the linear response is a combination of the responses in different natural modes. (In view 
of Fig. 14, for the base shear, T = 0.1 sec.) Consequently, we can expect the good performance of the SEB 
THAAT and Eq. (11) when applied to the response history analysis of the structure, especially for the top 
displacement. It is meanwhile worth noting that, different from the previous examples, where the damping was 
classical viscous (in the first example) and zero (in the second and third examples), in view of Fig. 13(a) and 
Table 10, the damping is non-zero and non-classical viscous in this example [2,36,37]. 

The model is analysed twice using the average acceleration method [33], and twice using the 
Wilson-@ [38,39] (@ = 1.4) method. In the first analyses by either method, the integration step is set to -At 
(obtained from Eq. (2) for the first target response). The analyses are then repeated using At = 2 ,At, after 
implementation of the SEB THAAT. The fractional time stepping method [40,41] is used for nonlinearity 
solution, and the nonlinearity tolerance and the maximum number of iterations are considered equal to 1E-6 
and 5, respectively (see [17,40,41]). The results are displayed in Fig. 16, and are evaluated, taking into account 
that the responses of nonlinear dynamic analyses may be inaccurate even significantly regardless of the SEB 
THAAT [1,3,7,28,42-46]. For the first target response, i.e. the top displacement, the responses obtained 
with or without application of the SEB THAAT coincide, when analysing with either integration method. 
For the base shear, the two responses are close, in analysis with the average acceleration method [33]. In 
analysis with the Wilson-6 [38,39] (@ = 1.4) method, however, the computed two base shears are evidently 
different. To better study the difference between the two base shears, obtained from_ the 
Wilson-8 [38,39] (@ = 1.4) method (see Fig. 16(b)), the results of analysis with very small steps is displayed 
in Fig.16(c). In view of this figure, both of the base shears computed by the Wilson-9 [38,39] (@ = 1.4) 
method (when applying and not applying the SEB THAAT) differ significantly from the exact base shear 
(displayed in Fig. 16(c)). In more detail, the two differences (between the two base shears in Fig. 16(b) and 
the exact base shear in Fig. 16(c)) are much larger than the difference between the two base shears in Fig. 
16(b). Therefore, the accuracy of the base shears in Fig. 16(b), is acceptable, in the sense that replacing the 
base shear obtained using At = ,At with the base shear obtained using At = 2 ,At (both displayed in Fig. 
16(b)) does not imply meaningful change in the response accuracy. Even more, returning to the origin of the 
Wilson-8 method [38], one of the main purposes of the Wilson-@ method is to filter high mode responses out 
of the response [47]. This filtering, which is broadly known as numerical damping of time integration methods 
[3,4,6,20,21,48,49], is essential, when we are seeking the responses of the structural model before the 
discretization resulting in Eq. (1). The discretization, replaces the structural model (with infinite number of 
degrees of freedom) with the mathematical model in Eq. (1) (with finite number of degrees of freedom), in the 
price of spurious high frequency oscillations in the response, which can be eliminated by numerical damping 
[3,6,49]. (Numerical damping can also eliminate real high modes with small contribution in the response, but 
computed erroneously because of largeness of the At/T [3,6,47-49]). Considering this, when using the 
Wilson-@ method [38], the purpose of the analysis may be different from achieving good accuracy compared 
to the exact response reported in Fig. 16(c). As a result, because of “selection” of the integration method, the 
base shears in Figs. 16(a) and 16(b) not only show the good performance of SEB THAAT when using Eq. 
(11), but can also be considered sufficiently accurate. 
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Fig. 16. Target responses of the fourth example: (a) Computed using the average 
acceleration method, (b) Computed using the Wilson-6 (@ =1.4) method, (c) Exact 
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For further clarity, the following two questions are to be answered, as well: 

(a) What is the reason of the different performances of the SEB THAAT using Eq. (11), in Figs. 16(a) and 
16(b)? Meanwhile, why are the base shears in Figs. 16(a) (and Fig. 16(c)) and 16(b) such different that 
while in Fig. 16(a) (and Fig. 16(c)) the base shear increases with time, in Fig. 16(b), the base shear 
decreases with time? 

(b) What is the reason of the negligible difference between the top displacements in Fig 16(b), while the 
difference between the two base shears in Fig. 16(b) is recognizable? 

The answer to Question (a) lies in the fact that Eq. (2) is not rigorous; besides other ambiguities, it is not 
clear why the integration method does not affect the selection of the integration step. More specifically, the 
numerical damping referred to in the previous discussion is very different in the average acceleration and 
Wilson-8 methods, with no influence on Eq. (2). The difference is evident in Fig. 17, for problems with 
classical viscous damping [50], and is yet unstudied for problems with non-classical viscous damping. (In Fig. 
17, p and T imply the spectral radius [6,20,48-50] and the oscillations period, respectively, and the difference 
of spectral radius from one at large value of At/T represents the capability of eliminating the oscillations with 
period T in analysis with step At.) Despite the latter, there are experiences in the literature (e.g. see [51]), for 
extending the discussions on classical viscous damping to non-classical viscous damping, using complex 
variables [52], such that the classical case can be considered as a special case of the whole discussion. 
Therefore, it is reasonable to consider significant difference between the numerical damping of average 
acceleration and Wilson-@ methods, when the viscous damping is non-classical, as well. This explains the 
considerable difference between the base shears in Figs. 16(a) and 16(b). The numerical details however cannot 
be presented, because the figure corresponding to Fig. 17 is yet unavailable for problems with non-classical 
viscous damping (the case in this example). Only, as a simple rough study, the mid-parts of the base shears in 
Figs. 16(b) and 16(c) are compared in Fig. 18. Accordingly, the change of the base shear (versus time), 
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Fig. 17. Changes of spectral radius (p) with respect to At/T for the 
(a) average acceleration method [33], (b) Wilson-6 (0 =1.4) [38,39] 
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Fig. 18. A simple rough comparison between the base shears in Figs. 16(b) and 16(c): 
(a) Middle ten seconds of Fig. 16(c), (b) Middle ten seconds of Fig. 16(b) (the left figure), 
(c) Middle ten seconds of Fig. 16(b) (the right figure) 
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including the change of base shear increase (in Fig. 16(c)) to the base shear decrease (in Fig. 16(b)), occurs 
along with further removal of high frequency oscillations, in analysis with larger integration steps. 

In answer to Question (b), if we compare the exact top displacement and the exact base shear in Fig. 16(c) 
(or even in Figs. 16(a) and 16(b)), the value of T is much larger for the top displacement. As a result, the value 
of At/T and accordingly the value of At/T is much smaller for the top displacement, compared to the base 
shear. From the other side of view, because of the essentiality of convergence for time integration methods 
[3,4,6,11,12,20,21], the difference between the spectral radii of different time integration methods disappears 
for sufficiently small values of At/T. (This is evident for problems with classical viscous damping in Fig. 17, 
as well as in the extended study considering many integration methods reported in [50].) Consequently, the 
effect of the numerical damping of Wilson-@ method on the top displacement is much less than the effect on 
the base shear. In other words, the integration has removed the high frequency oscillations from the response, 
and since the high frequency oscillations have negligible contribution to the top displacement (see Fig. 16(c)), 
the two top displacements in Fig. 16(b) are different negligibly. Due to a similar reason, the difference between 
the two base shears in Fig. 16(b) is noteworthy. This plus the role of T in Eq. (2) (though not rigorous) 
completes the answer to Question (b). 

Finally, the reductions in analysis run-time, due to the SEB THAAT using Eq. (11), are equal to 29.32% 
and 22.34%, when using the Wilson-9 (6 = 1.4) and average acceleration methods, respectively. The 
difference of these values with 50% is because of the nonlinearity of the problem, explained previously in this 


paper. 
Complementary discussion 


In the previous sections, it was demonstrated that, in seismic response history analysis of mid-rise steel- 
structure buildings, by using Eq. (11), the SEB THAAT can reduce the analysis run-time, leading to 
sufficiently accurate responses, without prior knowledge about the response. Besides, in practice, buildings 
are subjected not only to ground motions, but also to gravity loads. As a result, the effect of using the SEB 
THAAT and Eq. (11) on the real analyses is even better than that discussed in the previous sections. This 
enhances the importance of the simplicity obtained from Eq. (11). Further discussion on some ambiguities and 
limitations is however essential. 

The first ambiguity is whether the good performance observed in the presented examples relates to the 
seismic design code. In other words, is the observed good performance limited to the analysis of structures 
designed using the Iranian seismic standards? By using different design codes, the results of response history 
analysis with/without using the SEB THAAT will change. The simplicity and good performance of using the 
SEB THAAT and Eq. (11) will however not change, because of two main reasons. First, the scientific bases 
of seismic codes are close, and hence for a similar seismicity, the designs obtained from using different codes 
generally differ slightly. As a result, the effect of the change of the code on the performance of the SEB 
THAAT taking into account Eq. (11) would be reasonably tolerable. The second reason is that the codes used 
in the presented examples were the standards of Iran, and the seismicity of Iran is higher than many other 
regions of the world [53]. Therefore, when changing the design codes, the resulting T in Eqs. (2), (3), and (5) 
will probably increase. This implies even more suitability for Eq. (11) and more accuracy for the computed 
responses, when using seismic codes different from the codes of Iran. 

In view of the presented examples (see Figs. 5 and 11), the second ambiguity is whether the SEB THAAT’s 
performance is acceptable when the structural system is irregular in height or plan. In two separate studies 
[54,55], attention was paid to the performance of the SEB THAAT when applied to analysis of buildings’ 
structures with irregularity in height or plan. In both studies, Eq. (11) provided an appropriate selection for n. 
A similar observation was made in a slightly different study on the SEB THAAT [56], as well. Accordingly, 
applying the SEB THAAT to seismic response history analysis of mid-rise steel-structure buildings with 
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irregularity in height or plan may be successful, for values of n equal to or greater than two. Irregularities 
simultaneous in both height and plan are yet not tested for the performance of the SEB THAAT. Considering 
this and for the sake of brevity, none of the tests on irregular structures is reported here. Therefore, the claims 
in this paper are limited to mid-rise steel-structure buildings categorized regular by the seismic codes. This 
limitation is however practically unimportant, because, it is decades that buildings’ designers mostly prefer to 
design the structures to be regular according to the seismic codes!??0?!2:23.4 ; see also [25,26,57]. 

In view of the presented discussions, it is notable to add that few successful tests are reported on taller 
buildings [58] and buildings with concrete structure [56], as well. Besides, in view of the presented discussions 
and examples, and the more studied examples, not reported here for the sake of brevity, no limitation seems 
existing on the time integration method. Further investigation is essential. 

Considering the limitations addressed above, the social importance and large number of mid-rise buildings, 
the simplicity of Eq. (11), the everyday smaller values of digitization steps, the generally time-consuming 
nature of response history analyses, and the considerable reductions in run-time reported in the presented 
examples, the achievements are significant. Therefore, the future of the presented research is promising, at 
least until Eq. (11) can be replaced with a better comment or computational procedure (see also [10]). 


Conclusion 


The SEB THAAT is a technique for accelerating different analyses of structural systems. In this paper, an 
engineering comment for applying the SEB THAAT in seismic response history analysis of mid-rise steel- 
structure buildings without any details about the response is proposed. The engineering comment implies 
obtaining the step enlargement scale n, in application of the SEB THAAT to response history analysis of mid- 
rise steel-structure buildings, from Eq. (11). By implementation of this comment, without notable effects on 
the response accuracy: 

(a) The SEB THAAT can be applied to response history analysis of mid-rise steel-structure buildings, in a 
much simpler way (Step 6 in Fig. 2 is considerably simplified). 

(b) The analysis run-times can be significantly reduced (50% for linear analyses). 

This is a significant achievement, which, as a secondary achievement, can encourage structural analysts of 

mid-rise steel-structure buildings to use response history analysis in real projects. 

Limitations exist for the proposed comment. In addition to those implied in the expression “mid-rise steel- 
structure buildings”, the most important limitation is that the building structure must be “regular” according to 
the seismic code. This is however not a severe limitation, in view of the results of some recent researches, and 
the current practice of buildings structural design. Besides, the accuracy of the obtained responses need to be 
interpreted considering the numerical damping of the integration method. Other minor limitations exist, 
as well. 

The future of the proposed comment is promising, with attention to its simplicity and effectiveness, and the 
fact that digitization steps are in every day decrease. Finally, extending the application of the proposed 
comment from response history analysis of mid-rise steel-structure buildings to other classes of analyses, and 
improvement of this comment are two areas for further research. 


'S BHRC (Building and Housing Research Centre), Standard No. 2800-05, Iranian Code of Practice for Seismic Resistant 
Design, Iran, 2007 (in Persian). 

20 INBR (Iranian National Building Regulations), Iranian National Building Code, Part 10-Steel Structures, Iran, 1993 
(in Persian). 

21 ICC (International Code Council), IBC - International Building Code, Club Hills, USA, 2003. 

22 BCJ (Building Centre of Japan), Structural Provisions for Building Structures, Tokyo, Japan, 2001. 

23 RAK, Greek Code for Seismic Resistant Structures, Athens, Greece, 2000. 

24 NRCC (National Research Council Canada), National Building Code of Canada, Canada, 2005. 
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but using the current techniques only reduces the volume of water 

provided to the zone by 8 to 10 percent and the time of supply by a fixed 

amount (not around the clock). Additional investigation led to finding © The Author(s) 2022 
the causes of zoning's poor effectiveness, the shortcomings of the 
design, insufficient level of network research, application of incorrect 
principles of zoning and pressure management. We developed the 
methods discussed in this paper to resolve the aforementioned Tie spel ts licsieed made aCe 
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until the night-time consumption of the zone is less than half of the 

daytime consumption. 

Keywords: water supply network, water supply network zoning, 
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Introduction 


The requirement for zoning of the distribution network arises in mountainous conditions when there is a 
significant difference in levels within the boundaries of the settlement. Both gravity and pressing systems 
(reverse zoning) may require zoning [1,2]. 

When high-rise multi-residential buildings are constructed in areas where low-rise buildings are already 
present or when the high-rise buildings’ absolute height exceeds that of the low-rise buildings within the water 
supply zone's boundaries, zoning becomes crucial. Therefore, to increase the controllability of the distribution 
network, it is necessary to implement effective zoning of the network, i.e., to transform the network into 
hydraulically separated zones that are isolated from one another. Isolation must be implemented by installing 
existing or new valves, which will be closed during the regular operation of the water supply network, but can 
be opened if necessary. In mountainous areas, it is advisable to use vertical zoning in a sequential or parallel 
scheme. 

Many automatic network zoning approaches have recently become widely used with specifically created 
computer programs (hybrid approaches for the automatic partitioning of a water distribution network, spatial 
analysis zoning approach, semi-supervised method etc.). However, the application of these techniques is 
unworkable with an existing network that is worn out, elementally constructed and reconstructed, or 
fragmented. Our original research showed that some additional problems arise when creating zones using 
existing weathered, prefabricated, reconstructed, or fragmented networks. Thus, after zoning Arabkir zone 1, 
Kievan, Gulbekyan, and Keri streets, it was found that the amount of water supplied to the zone decreased by 
only 8%—10%, and the 8-hour duration of supply was 12 h instead of the planned round-the-clock. 
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Because the designed pressure wasn't maintained, accidents happened more frequently. Additional 
investigation led to finding the causes of zoning's poor effectiveness: 
e the shortcomings of the design, 
e insufficient level of network research, 
e application of incorrect principles of zoning and pressure management. 

We developed the methods discussed in this paper to resolve the aforementioned problems. 

Because of the analysis of the original research and operation data on the water supply network, as well as 
the international experience applied in similar conditions and the study of the existing literature, the zoning 
works carried out in the water supply network of different settlements of RA have been improved and are 
currently being carried out efficiently according to a procedure developed by us. The principles of water supply 
network zoning are found in several versions of the literature currently in print. Many automatic network 
zoning approaches have recently become widely used with specifically created computer programs. However, 
the application of these techniques is unworkable with an existing network that is worn out, elementally 
constructed and reconstructed, or fragmented. 

Different models have been proposed, based on classic optimization methodologies or on meta-heuristic 
approaches [3-7]. 

The optimal design of DMAs as part of a Decision Support System (DSS) for reducing the water losses has 
been addressed only recently in the literature. Some authors have proposed hybrid approaches for the automatic 
partitioning of a water distribution network, based on both meta-heuristic algorithms and on applications from 
the graph theory [8,9]. 

The proposed models mainly focus on the preservation of the hydraulic reliability of the network, while 
less control is allowed on the costs of the provided solutions [10]. 

Sempewo et al. [11] developed a spatial analysis zoning approach based on the METIS graph partitioning 
tool [12]. 

The proposed technique follows the analogy with the distributed computing methodology of equally 
distributing workloads among processors. However, as stated by the same authors, although the method results 
effective in the demarcation of contiguous districts, the quality of the provided solution degrades when 
considering multi-objective partitioning, and uncertainties are produced when increasing the number of DMAs. 

Herrera et al. [13] proposed a semi-supervised method named multi-agent adaptive boost clustering. This 
complex technique considers both the WDN features (e.g. node elevations and demands) and the economic 
issues (edge cut = number of pipes to be intercepted) for the partitioning of the network. Nevertheless, the 
procedure is only applied to cases in which the number of DMAs is lower than the number of supply nodes. 

More recently, [14] have introduced an approach for automatic creation of DMAs based on the hierarchical 
community structure of the WDN. In this study, the selection of the feeding lines for the DMAs is made through 
an iterative selection method based on a sensitivity analysis. The methodology was tested on a very large 
network with reasonable computing times, but the results showed high sensitivity towards the assignment of 
the input parameters. 

A comprehensive description of the possible objective functions for the problem can be found in 
Gomes et al. [15]. 

The two-step approach proposed by the authors consists of a preliminary partitioning of the WDN into 
suitable DMAs through the application of the design criteria and graph theory concepts, i.e. the Floyd-Warshall 
Algorithm (FWA). Next, Simulated Annealing (SA) is used for the localization of entry points and of boundary 
valves, but also for planning the reinforcement/replacement of pipes in the network. Although the results have 
resulted satisfactory, the global optimality of the solution is not ensured by the SA. Hence, despite the newest 
contributions, further developments seem required about this topic [16]. 
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Among the available approaches for pressure management, the use of District Metered Areas (DMAs) also 
allows for a more accurate localization of the leakages in the water distribution network, which is achieved by 
monitoring the input and the output discharges for each district. Nowadays, this approach is widely used in 
practice, but its application is still largely entrusted to the experience of technicians [10]. 

According to the proposed principles, to increase the zoning and manageability of the water supply network, 
it is necessary to create hydraulically isolated zones with up to 3000 subscribers in the settlement. Since the 
zoning is carried out under the conditions of the existing system, the procedures developed during the design 
and reconstruction work consider the features of the existing system: difference in characters; dictating 
characters of the day regulatory reservoirs; location and volumes; diameters of the existing pipes; height of the 
buildings; population density; etc. 

It is suggested to remove the inner rings of the zone by temporarily converting the ring network within the 
presented zone boundaries into a dead-end network using valves and separate sectors. Magnetic flowmeters 
and pressure recording sensors should be installed on the feeding pipes. It will allow for assessing the water 
balance within individual zones, determining the amount of unaccounted water, and discovering its nature and 
location [17]. 

These investigations aim to calculate how much water enters each zone and compare it to how much water 
is received based on the number of subscribers. By creating water supply sectors with one-way flow and an 
exact number of subscribers, it is possible to estimate the amount of unaccounted water in the presented zone. 
During the zoning (design and construction) of the existing water supply network, it may be necessary to apply 
the following procedures: changes of water sources (spring, aqueduct, daily regulating reservoir) in zones or 
sometimes in entire districts, change of design and/or operating hydraulic regimes of water pipelines, 
converting the ring network temporarily into dead-end sectors with separate one-way flow, loss detection and 
elimination, installation of pressure regulators at required points, decommissioning or replacing backyard 
pumping stations with smaller capacity pumps, as well as a grouping. 

Before the implementation, it is crucial to perform studies and designs at a high level through experienced 
engineering solutions because zonings require financial investments and impact the quality of services offered. 
However, experience shows that it is almost impossible to have a final zoning project, so during the execution 
of the works, there is a need to review the scope of the study. 


Materials and Methods 
Design and implementation of a water supply zone 


Before starting the zoning, it is necessary to analyze the problems and data recorded during the operation, 
define the registration procedure, create a working group, acquire the needed equipment for measuring and 
control, and eliminate obvious (visible) leakages. 

The following sequence of work has been developed for the design of zones and their implementation 
(construction). 

1. The existing water supply distribution network is divided into separate zones through valves on the plan. 
During this demarcation, it is necessary to consider the possibility of creating optimal pressure regimes in 
the zones and excluding the dead sections in pipelines. 

Limiting valves should be installed just away following the connections of relatively large volume of water 
users to prevent the creation of dead-end areas (no flow or small flows), as the valves turn the pressure 
zone into a dead-end network. Water quality degradation is possible in areas with a low water flow rate 
[1]. As mentioned above, to ensure the controllability in the pressure zone, it is necessary to create sectors 
with a smaller number of subscribers. According to the results of the studies, the zone should include 
500-3000 subscribers, depending on the level of leakage, the development characteristics of the district, 
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the applied methods of leakage control, and hydraulic conditions. The water supply zone may serve 3000 
customers in densely populated places, as it is in the settlement’s city center. However, in these cases, it is 
challenging to distinguish minor breaks separated from flow data recorded at night, making it more difficult 
to pinpoint their place. However, by temporarily closing valves, large water supply zones can be divided 
into sub-sectors of two or smaller sizes. In this situation, extra valves might need to be placed during the 
water supply zones’ design phase. 

2. The power source is determined based on the case of providing an uninterrupted water supply to the facility 
in the most dangerous (critical) conditions of the specified area. In design practice, these zone feeding 
options are possible: 

1) the zone may have one or two power sources: an aquifer or a reservoir, 
2) zoning work can be carried out sequentially: here, the second zone receives water from the first zone 
in a transit way (Fig.). 


2" pressire zone 


1“pressire zone 


daily regulatory 
reservoir 


Fig. Sequential zoning scheme 


3. It is planned to close off the valves on all border feeding lines during the hydraulic isolation of the divided 
zones, leaving only the targeted feeding line. To guarantee total isolation, sometimes required to install 
new valves. Small-scale maps of the distribution of water lines, information from operating staff, and 
existing hydraulic data are employed at this stage of contour limitation of the water supply zone design. 
In establishing the boundaries of the water supply zones, besides the general design criteria, the following 
conditions should be considered: 

e make the most of the feed source's power or pressure (day regulating reservoir, aqueduct) position. In 
this case, area and topography may dictate sequential zoning using existing pipelines. 

e any transit aqueducts or water lines feeding other zones are avoided as much as possible while drawing 
the boundaries of a zone, 

e if possible, the diameter of the water line supplying the zone should not exceed 300 mm to avoid high 
financial costs and technical complications associated with the flow measurement unit installation, 

e to ensure that research can take place, the area should have a constant supply of water. 

4. "Zero tests" is used to verify the compliance of the design works with the completed reconstruction works. 

It provides the implementation of the following conditions: 

e inspecting the installed flow meter for accuracy, 

e checking the functioning (hermetic) of the boundary valves, 

e detection of hydraulic connections between neighbouring zones, 

e inspection and assessment of changes in the quality of water supplied to customers in the reconstructed 
zone and neighboring zones, 


e check for excess pressures. 
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Zero pressure test procedure 


Zero pressure testing is carried out at night, between 1° and 5°. In order to detect hidden connections 
between zones, data recording devices (loggers) are connected to the cost meters of the lines feeding the tested 
and neighboring zones. 

The tightness of the boundary valves and the existence of any hidden connections are then assessed by 
closing the boundary valves. For the situation assessment, it is also necessary to temporarily install the 
automatic pressure recording devices at defined points in the network during testing. After that, the valve 
supplying the zone is closed, and if there is no significant pressure drop within 5 minutes, any hydrant in the 
lower part of the network is opened, reducing the pressure to zero. The pressure in the water supply zone must 
be zero when the hydrant is closed; otherwise, another power source that wasn't considered during the design 
may exist, or the tightness of any boundary valve that is audibly checked may not be maintained. 

The zero pressure tests are a crucial procedure for determining whether the water supply zone is 
hermetically sealed, so the aforementioned investigations must be continued until the experiment’s successful 
completion is guaranteed. It may take up to a year to get positive results from the zero pressure test, depending 
on network characteristics (zone size, number of subscribers, population density, and the ratio of private and 
multi-residential buildings). 

The described measures are necessary because proper zoning was not carried out during the design and 
construction of existing water supply networks. Besides, the circumstance that the drawings created during the 
construction of the water supply system do not match the actual situation made the design and implementation 
of the zone more challenging. 

Failure detection and elimination are carried out to reduce power losses. For this phase, we determined that 
the night-time costs of the zone should not exceed 50% of the daytime cost. Finding and eliminating hidden 
leaks can also take several months. 

Day and night flow is recorded after these functions' completion and is necessary for benchmark indicators 
to evaluate the effectiveness of zoning (redevelopment), which must be monitored and maintained during 
operation. 

The precise estimation of the zone's night flow and its ongoing regulation, as well as data updates, are of 
particular importance required for the following functions: 

e continuous maintenance of the water supply zone's hydraulic isolation. There is often a conflict of interests 
between the operating units and the zoning group: opening of boundary valves, disturbance of zone 
boundaries, etc. Such cases must be monitored and prevented until their elimination, 

e regularly updating the number of subscribers, 

e modernization and improvement of data acquisition methods and technologies, 

e leak detection and elimination in a planned manner, 

e regular analysis of received data. 


Conditions for creating a sector (dead - end network) during the design 


Because of the formation of sectors within the borders of a separate zone, the existing ring network turns 
into a dead-end having a one-way supply. Hydraulic calculations are required to verify the transmittivity and 
pressure losses of the sections to maintain the water supply of the dead-end network. 

For this purpose, a calculation scheme for the dead-end network has been created. It is divided into 
calculation sections, where the directions of water flow are noted, and the actual outputs in the sections are 
measured. 

In contrast to the well-known method for calculating the dead-end network, during the conducted studies, 
the outputs obtained through experimental measurements are taken as the calculation output for each section 
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based on the maximum actual water demand. When choosing the endpoint for constructing pressure lines based 
on the measured outputs, not only the distance and relief but also the required pressure at that point, considering 
the height and location of the buildings, are considered. 

After measuring the outputs of the designed dead-end network segments, the network transmittivity 
condition is checked. For this purpose, the pressure losses occurring in the sections are determined by the 
formula h)= SQ’, using F.A. from Shevelyov's tables [3], based on which the free pressure lines of the network 
are built to have the pressure magnitude at all points of the selected calculation direction. Free pressure lines 
are installed in every direction possible to assess the pressure at the examined network's sites. In the event of 
a problem providing pressure in the subzone of the designed dead-end network, based on specific conditions, 
it can be solved by the following options: 

e by over-correcting the pressure regulators or opening the compressed valves on the supplying water pipe 
(if available), 

e the excess pressure regulation increased the pressure in the initial part of the considered sub-zone in the 
neighboring sub-sector, 

e to change the power source of the zone that will provide the required pressure or apply another zone 
scheme, 

e pressure losses in the specified sections reveal the significant leaks by finding and fixing leaks in any 
section or sections. Therefore, it may be found that the cause of the pressure loss is a hidden local 
resistance (unknown compressed valve, presence of a gasket, or blockage), 

e by increasing the diameter of a specific network section if the pressure line in that section has a steep 
slope, 


e by installing local pumps if the pressure is insufficient for a few high-rise buildings. 


The last two options are advisable to use in case of economic feasibility. It should be added that when 
performing zoning, the need to increase the diameter or install a pump arises in rare cases because the diameters 
of pipelines built in the Soviet era are chosen with a high stockage. 


The priority of the implementation of distribution network zoning 


Because of the zoning processes’ financial and technical requirements, it is impossible to carry out 
operations simultaneously in all planned zones. Instead, priority is given to the zones with the highest leakage 
levels. Based on the technical condition of the distribution network of the settlement, it is recommended to use 
the following expression of specific losses of flow: 


qleak = (Qnight = Ocous) / Lw ( 1/hour m), (1) 


where 
leak - Specific loss of flow, I/s m, 
Qnight - average night flow recorded by the zone flowmeter, I/h 


Qcons, - the average consumption of subscribers during the night I/h, 


Ly - length of pipelines of the zone, m. 


In the sector of own residences, the accepted average consumption per subscriber during the night was 1.7 
I/h, while in multi-apartment buildings, it was 0.6 l/h. These numbers are used to calculate the amount of 
leakage in the systems of European cities where operational circumstances match up with regional water 
supply systems. 

Because customer overnight costs can significantly affect the zone's estimated night loss, they are calculated 
independently using actual data. 
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Results and Discussion 


Analyzing the results of multi-year studies and considering the current technical situation of the distribution 
network, it is recommended to continue the accident detection and elimination work for the current phase until 
the night-time consumption of the zone is less than half of the daytime consumption: 

Qnight = Qaay /2 (1/s) > (2) 
where Qnight and Qqay are the average water volumes given to the zone during the night (1°°-5° period) and 
daytime. 

The mentioned expression was defined after studying the zones with the lowest water losses in RA 
settlements (Davtashen, South West, and other districts). Thus, let's discuss the study results on the 
"Arabkir 1" zone, where the average amount of water supplied at night is 102 I/s, and during the day is 133 I/s. 
Applying the expression (2), we can determine the amount of leakage reduction in the “Arabkir 1” zone to 
meet the current stage requirements, which should not be greater than 133/2=66.5 I/s. From the data obtained, 
we can conclude that in the mentioned zone, there is still leakage of 102-66.5=35.51/s more than the permissible 


one, which needs to be detected and eliminated. 


Conclusion 


To increase the efficiency of the aqueduct network, zoning problems were studied, taking into account the 
technical condition of the system. Scientific calculations were carried out, and a methodology for zoning was 


developed: 

e The analyses showed that the water pipe network should be transformed into hydraulically separated 
zones while using the current distribution system, taking into consideration the position, placement, and 
dependence of the site's relief features. It was necessary to create a zone design and implementation that 
considers the principles of system management, the formation of sub-zones (dead-end networks), the 
calculation of optimal pressures, the detection and elimination of current losses, etc. 

e To increase the level of management of zone distribution networks, the means of automatically 
transferring the fundamental characteristics of operation: pressure, output, and electricity consumption 
values are applied, the investment of which in water supply systems will give the desired results. 
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Introduction 


Operation practice has shown [1,2] that two automobiles of the same brand have different levels of a 
residual resource after a particular mileage’, both in terms of constituent units and the whole automobile [2]. 

It means that the automobile physical wear coefficient [3] needs to be considered more 
deeply by analyzing and taking into account some technical and technological factors during the operation, 
which forms the rolling stock life cycle in the given operating conditions [4]. The 
proposed concept for determining the automobile depreciation period is based on maintaining the smooth 
operation of the rolling stock throughout its entire life cycle. Technical disruptions that cause the transport 
process to stop, result in inefficient downtime as well as material and labour expenses [1]. To avoid accidental 
and sudden disruptions, the actual coefficient of the automobile's physical wear? should characterize the current 
technical condition of the rolling stock and engine constituents, ensuring the need for repair and maintenance 
within a given period of millage and time [5,6]. Up to date, the current methods for determining the operating 
life have been of analytical nature. They do not reveal how the physical wear coefficient of the automobile 
varies due to the increase in overall mileage, exact operating conditions, and other factors [4,10]. 

The paper proves the stochastic nature of variations in the automobile's physical wear 
coefficient, it is considered as a random value. The characteristics of its variations pattern have been 
determined, allowing the automobile depreciation period to be calculated based on_ the 
economic indicator value preferred by the given economic entity. 


Materials and Methods 
It is necessary to develop qualitative and quantitative indicators for assessing the automobile’s physical 


wear aimed at solving the proposed problem and then discovering the analytical connections of their 


' GOST R 50779.10-2000 
> R 03112194-0376-98. Metodika otsenki ostatochnoy stoimosti transportnykh sredstv s uchetom tekhnicheskogo 
sostoyaniya, 1998 (in Russian). 
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interactions and relations through theoretical and scientific-experimental research [6]. It is well-known that the 
resource indicators of automobile units, auto parts, and interchanges have quite different values and are of 
stochastic nature [7]. A scientific-experimental research has been conducted for a group of minibuses (30 units) 
performing intra-city passenger transportation to study the variation in the automobile physical wear 
coefficient and identify the distribution patterns as a random value’. 

The following indicators have been studied during the scientific-experimental research: 

The number of automobile disruptions (nd) by years and hence the number of downtime (ndtp) 

days caused by them, 

The average daily mileage of one automobile (fd) —per km, 

The average annual mileage of one automobile (fa), per thousand km, 

The total mileage of the automobile group (>) £) per thousand km, 

The number of annual operation days of one automobile (nw) per day, 

The coefficient (a) of technical readiness of the automobile group, 

The disruption flow parameter (@) of the automobile group per-thousand km, 

The average mileage of automobile smooth operation (£) per thousand km fs, 

The physical wear coefficient (K) of the automobile group. 

The graphs (Figs. 1, 2) based on the data presented in the Table indicate the variations in the coefficient of 
the automobile technical readiness, average mileage of smooth operation, number of maintenance downtimes, 
and physical wear coefficient. According to the presented data, the coefficient of the automobile’s technical 
readiness during the fifth year of operation was 0.73, which is less than the average value of this coefficient 
(0.8). A similar solution is observed in terms of the average mileage of the smooth operation of automobiles. 
In case of an average value of 12.8 thousand km, the indicator for the fifth year of operation was 8.87 
thousand km or decreased by 3.93 thousand km (30.9%). 


Fig. 1. Dynamics of quantitative and qualitative variation of automobile 
operating indicators (a (1), €s. (2), ndt. (3)) by years 


3 GOST R 50779.10-2000 
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Fig. 2. Change in the automobile physical wear coefficient (K) and the average number 
of working days per automobile (nw) according to the years of operation 

Further analysis of the automobile’s smooth operation shows that, according to the operation results of the 
sixth-year, the indicator was 3.03 thousand km. The same dynamics is observed in the operation results seven 
year. According to the data, the automobile’s technical readiness coefficient during the fifth year of operation 
was 0.73, which is less than the average value of this coefficient (0.8). A similar situation is with the smooth 
operation of automobiles. In the case of an average value of 12.8 thousand km, the index for the fifth year of 
operation was 8.87 thousand km or decreased by 3.93 thousand km (30.9%). It indicates the following: if we 
are to be directed by the dynamics of the change in the average mileage of automobiles, we must decomission 
the rolling stock at the end of the fifth year of operation or the beginning of the sixth year due to low technical 
and economic indicators. 

Table. Operational indicators of the Automobile group 


N Year of operation 
1 2 3 4 5 6 7 
Indicator 
Number of disruptions na 60 78 92 110 148 388 516 
2. | Downtime na (automobile-per day) 209 750 1134 1726 | 2894 3557 4804 
TS — 2 (number) 203 164 140 127 101 91 64 
Current repair (number) 6 586 998 1599 2793 3466 4740 
3. | Working days, annual average(nw) 358 340 327 307 268 246 204 
Total 10741 | 10200 | 9816 | 9224 | 8056 | 7399 6146 
4. | Coefficient of technical readiness, a 0.98 0.93 0.89 0.84 0.73 0.67 0.56 
5. | Average daily mileage, fa (km) 218 210 186 179 163 160 135 
6. | Total mileage, } £ (thousand km) 2341.5} 2142.0} 1825.8} 1641.1} 1313 | 1182.9} 829.7 
7. | Average mileage of smooth operation, | 44.5 27.3 19.85 15.0 8.87 3.05 1.61 
£, (thousand km) 
8. | Disruption flow parameter 0.0225) 0.036 | 0.050 | 0.0666} 0.113 | 0.328 ) 0.578 
w(disruption per thousand km) 
9. | Physical wear coefficient, K 0.023 | 0.039 | 0.056) 0.0748} 0.154 | 0.490 | 1.032 


34 


Y. Vardanyan, V. Harutyunyan, V. Koichev, K. Mosikyan 


Now the problem should be considered in light of the dynamics of changing the indicator of inefficient 
downtime of the automobile. 

The average indicator of inefficient downtime for a group of automobiles' downtime was 2150 automobiles 
per day, or 71.67 automobiles per day for one automobile. Moreover, during the fifth year of operation (Table, 
Fig. 1), the average downtime was 96.5 automobiles per day, which is about 25 automobile per day more than 
the average value. If we observe the variations in the number of automobiles’ operation days over the years, it 
becomes evident that in the fifth year, there are already 8056 automobiles per day, with an average value of 
8799 automobile per day. It equates to 35.9% inefficient automotive downtime in regard to the number of 
operation days. Naturally, the monthly number of automobiles operating by the order was 74.1%, which is 
inefficient from an operational standpoint (the accepted efficiency amounts to no less than 80% from a 
financial and economic standpoint). 


Discussion of Results 


To assess the efficiency of automobile operation, the physical wear coefficient should be observed on terms 
of quantitative and qualitative standards, as it is defined by inefficient dowtime spent on recovery and an 
indicator of recovery duration [3]. This is evidenced by the fact that the calculated physical wear coefficient 
according to the value of the automobile coefficient of technical readiness amounted to 0.154 in the fifth year, 
while in the sixth year it grew rapidly, amounting to 0.490. It means (Fig. 2) that the amount of inefficient 
downtime hasn't increased sharply, but the disruption flow parameter has increased, which has led to a sharp 
increase in the physical wear coefficient. If we consider the physical wear coefficient with the average millage 
of inefficient downtime and smooth operation, it turns out that this indicator has not undergone drastic changes. 

Based on the above-mentioned analysis results, it can be firmly insisted that the automobile's physical wear 
coefficient should be assessed from quantitative and qualitative perspectives, firstly - taking into account the 
quantitative value of the disruption flow parameter* and seondly - the amount of ineffective downtime days 
spent on the recovery from disruptions. 

To this end, the automobile physical wear coefficient should be considered as a random function, 1.e. 
according to the millage or years of operation of the rolling stock [6,11]. 

The automobile physical wear coefficient can be represented as a function of the rolling stock mileage f, 
in that case [6]: 

K =n(2). (1) 

At the same time, the automobile mileage is a variable with a probability density function of f (€), which 
causes variation of K parameters. Let use the Laplace transformation (1), which for the normal probability 
distribution law of random variables is [2]: 


f(@) = exp - eo |. 2) 


where o is the root-mean-square deviation, km (dispersion), tt - constant (3,14), @ - milage, km, 


7 ie milage mean path, km, and has the following form: 


p(z) = exp [-e ‘Zt a ; (3) 


2 
where Z is the constant function for normal distribution, @(z) is the function of K value. 

It can be assumed that the first two members of the function are the automobile physical wear coefficient 
K, and it is a linear function: 


K=a,ta,'f, (4) 


4 GOST R 50779.10-2000 
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where dg is the resource of the new automobile, a,- is the specific value of automobile resource change per 


mileage (1000 km). 
In this case, the function Laplace transformation (4) for the linear function is as follows: 
Px(Z) = Pe(a,*Z) exp(—ap *Z), (5) 
or 
ry 2.72 7472 
Qx(Z) = {exp |-?: ay'zt oa} exp[—a ‘zl, (6) 
Thus: 
= o-a,)?- Z? 
Ox 'Z = exp {l1- @ a1 +ag)e]+ "5 , (7) 


Comparing the expressions (7) and (3), it follows that the expression (7) is also considered a Laplace 
transformation for a new normal distribution value K, the average value of which is: 
KR=G540 <4. (8) 
where K is the mean physical wear coefficient. 
Root-mean square deviation is: 


Ox =Q,0, (9) 
and the coefficient of variation is 
1 
Ve = aT (10) 
ajo Vy 


The coefficient variation for the new random variable is dependent on the coefficient of variation of 
argument and the intensity of the parameter change, particularly on its increase. 

Considering the change in the automobile’s physical wear coefficient along with the increase in the total 
mileage (1), it becomes evident that the dynamics of the change in the physical wear coefficient has increasing 
nature. It means the probability density function of the automobile physical wear coefficient f (K), taking into 
account the expression (6) will have the following form: 


1 
f(kK)= aoVane*P [1 


(1) 


= A 
20?-a? 

It is also evident from the graph presented in Fig. 2, which comes to prove along with the automobile's 
physical wear coefficient (K), the number of operation days decreases as the maintenance downtime for the 
rolling stock repair increases. While solving the practical problems, if the function n(#) is not linear, but is 
continuous in the range of €; — @2 millage and is close to being linear, it can be replaced by a linear function 
[6]. 

Summarizing the results of the theoretical and experimental research, it becomes evident that the change in 
the physical wear coefficient under the specific conditions of an automobile operation has a stochastic nature 
with a normal distribution pattern. Considering the change in physical wear coefficient as a random value and 
determining the characteristics of its change pattern - mathematical expectation, dispersion, and variation 
coefficient, allow to determine the depreciation period of the automobile according to the value of the technical 
and economic indicator preferred by the economic entity. 


Conclusion 


The research on the changes in the automobile’s physical wear coefficient in the actual operating condition 
allows to identify the pattern of actual change of the coefficient, its mathematical expectation, root-mean- 
square deviation, and variation coefficient, as well as allows to assess the current value of the indicator in a 
definitive range. It will be possible if the density of the function f (K) is calculated, which in fact is related to 
intensity of operation (ag and a,), and it is the main and dominant factor in the change of the automobile’s 
physical wear coefficient, and the milestone of the developed methodology. 
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Abstract: The first U.S. caving expedition to Armenia, the South Caucasus, by Samvel Shahinyan* 
NSS cavers took place in August 2007. Subsequent expeditions took place in 
2010, 2011, and 2013, with additional trips planned for the future. The goal is 
the exploration and photo-documentation of the caves of Armenia, and to 


E-mail: sshahinyan @ nuaca.am 


increase awareness of its underground realms. Although, in the past, there had Received: 14.02.2023 
been a few known caving expeditions to Armenia, overall little information . 
existed. In addition, Armenia’s local caving community is small in number. As Revised: 03.03.2023 


a result, this topic was studied and the first official US caving expedition to Accepted: 18.03.2023 
Armenia was organized in 2007. During the first expedition, four of Armenia’s 

significant natural caves were explored in the province of Vayots Dzor: Mozrov 

Cave, Arjeri Cave (Cave of the Bears), Mageli Cave, and Karmir Cave (Red © The Author(s) 2023 
Cave). Man-made caves were also visited. Subsequent trips to Armenia in 
2010, 2011, and 2013 included (1) further exploration of Mozrov, Arjeri, and 
Mageli caves, (2) a cave trip to the neighboring independent Armenian 
Republic of Nagorno-Karabagh to explore Azokh Cave, and (3) the exploration ; as . 
of several caves in the northeast Armenian provinces of Tavush and Lori. This work is licensed under a Creative 
Natural caves consisting of limestone, conglomerate, and lava were explored See ae pane eae emeeal 
during these expeditions. Also, a number of man-made caves were visited, oO Bier ona TAcense 

some of which were used as churches in centuries past. This article summarizes 

the four expeditions and discusses both the natural and man-made caves of 

Armenia. I beleive the article will be interesting to builders, gas pipelines and 

road engineers. In practice, it can be used by travel agencies and individual 

tourists, as well as by all lovers of underground monuments of nature and 

culture. 

Keywords: expedition, survey, caves, cave plan, church, tourism, photo- 

documentation, cave map. 


Introduction 


The present-day Republic of Armenia lies geographically in the South Caucasus, between the Black Sea 
and the Caspian Sea. The most ancient of the countries and nationalities of the Caucasus, the indigenous 
Armenian people have a legacy stretching back 3.000 years. Eastern Turkey, southern Georgia, western 
Azerbaijan, and a northern portion of Iran are all part of historical Armenia. The Republic of Armenia is a 
rugged, mountainous land, with an average elevation of approximately 1.500 meters, and a population of just 
over three million (due to Armenia’s tumultuous and tragic history, there are presently more Armenians 
residing outside Armenia — the Diaspora). Armenia holds the distinction of being the first nation to declare 
Christianity as its state religion in 301 A.D. Ancient churches in Armenia (among the oldest in the world), 
many still standing and in use, pre-date those of Europe by centuries. Some of these churches were hewn out 
of rock, thus cave churches. 

Prior to 2007, available information and literature on natural and man-made caves in Armenia was sparse. 
After researching this topic and contacting local Armenians in Armenia, first United States (NSS) caving 


' A previous version of this paper was presented at the 17" International Congress of Speleology (Speleo 2017) 
in Sydney, Australia, on July 23 — 29, 2017. 
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expedition to Armenia was organized. The 
expedition was a success, as we explored and 
photo-documented four of Armenia’s 
significant caves in the province of Vayots Dzor 
in south central Armenia. Participating in this 
first expedition were Steven Johnson, James 
Wilson, Greg Chavdarian, Seda Chavdarian, 
and Charles Chavdarian of the United States, 
and Vrezh Nazarian and Samvel Shahinyan of 
Armenia, who also acted as our guides (Fig.1). 

This was followed by expeditions in 2010, 
2011, and 2013. Asa result, further exploration 
occurred not only in Vayots Dzor province, but 


also in the northern provinces of Tavush and 
ee ieee the eastern province of myunik: » Fig. 1. Members of the first U.S. caving expedition in 2007. 
addition, we also traveled to the neighboring James Wilson, Greg Chavdarian, Charles Chavdarian, Seda 
independent Armenian republic of Nagorno- Chavdarian, and Steven Johnson (Left to Right) — Republic of 
Karabagh. Both natural and man-made caves Ee ERO RY Vs MacNN 
have been explored and are discussed in this paper. The following cavers participated in some or all of these 
subsequent trips: Lara Chavdarian and Charles Chavdarian of the United States, and Vrezh Nazaryan, Smbat 


Davtyan, Pegor Papazian, and Nyree Abrahamian of Armenia. 


Materials and Methods 


Vayots Dzor Province, Republic of Armenia 


For the 2007 expedition, we maintained a base camp near the village of Mozrov, at an elevation of about 
1.700 meters. For subsequent expeditions in Vayots Dzor province in 2010 and 2011, we stayed in bed and 
breakfast homes in the town of Yeghegnadzor, the provincial capital. We explored and photo-documented 
several caves. 

Mozrov Cave: Mozrov Cave, primarily a limestone cave, was first discovered about 40 years ago [1] 
during road construction, when a collapse occurred resulting in the creation of a large entrance to the cave 
(which is also the only known entrance to the cave). The cave sits at an elevation of approximately 1.550 
meters adjacent to a mountain road. Due to its accessible location, the cave is vulnerable to visitation by non- 
cavers and tourists and has sustained some damage. Over the course of the four expeditions, we have explored 
the cave five times, and have taken a number of photographs. There is substantial decoration in this cave — 
stalactites, stalagmites, columns, moon milk, flowstone, coral, popcorn, crystalline spars, soda straws, 
helictites, etc. The cave decoration is also noteworthy for its myriad of colors — red, caramel, yellow, and 
white. We were to soon find that Mozrov Cave is not unique in this regard. The cave has over 300 meters of 
known passage. The cave consists primarily of an undulating large and long main chamber, with a separate 
large chamber that can be entered through a low and difficult to find passage at the back of the main chamber 
[2]. This second chamber is noteworthy for its extensive multi-colored decoration (and has been aptly named 
as “Decoration Hall”). 

Due to the ongoing damage that has been occurring within the cave, Mozrov Cave was recommended for 
tourism. As such, the cave would be gated and protected, and only controlled access would be allowed. Also, 
conservation of the cave would be instilled in the visitors. A management plan was created and provided to 
the foundation. To date, there has been no further action on this cave. It is our hope that Mozrov Cave will 
eventually be protected [3]. In 2010, at our request, our local Armenian caving colleagues, headed by Smbat 
Davtyan, returned and surveyed the cave and provided the first cave map of it (Figs. 2,3). 
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Fig. 2. Mozrov Cave Map — Survey 2010 — Smbat Davtyan 
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Arjeri Cave (Cave of the Bears): Arjeri 
Cave, at an elevation of nearly 1.700 meters, 
is Armenia’s largest known cave, with 
currently about 4 kilometers (2.5 miles) of 
passage. Our only map was a rudimentary 
overview of the cave which was created by 
Russian cavers nearly 30 years ago [4]. The 
cave does require an updated and detailed 
cave survey and map. A total of four trips 
were made into this cave in 2007 and 2011. 
Arjeri Cave, a limestone cave, is the most 


highly decorated cave in Armenia, and one of 


the most highly decorated caves we have ever Fig. 3. Mozrov Cave — Smbat Davtyan and Nyree Abrahamian in 
seen. To reach the only known cave Decoration Hall - Republic of Armenia. Photo by C. Chavdarian 
entrance, one must leave their 4WD vehicle 
off-road, and at an elevation of 1.500 meters, 
and undertake the final steep hike to the cave. 
Just inside the Entrance Hall, explorers are 
immediately greeted with the sight of huge 
flowstone columns (Figs. 4,5,6).After 
passing through this first large chamber, one 
enters a slope requiring a steep belly crawl 
upward over slick flowstone. Along the way, 
this area is decorated with various 
formations. This leads into a large chamber 
known as Photographer’s Hall, and the name 
is certainly appropriate as the room is rife 
with colorful formations —_ stalactites, 
stalagmites, columns, draperies, bacon, coral, 
and popcorn. The sheer quantity and density 
of the decoration in this large chamber is 
overwhelming. From here one continues on 
into Bear’s Hall, which contains the bones of 
a bear [2]. This, in fact, is the reason for the 
name of the cave. On the way to Vayk Hall, 
there is a deep pit, which is estimated to be at 
least 18 meters deep. This pit has not been 


explored by us (and may possibly lead to 
more passage). On entering Vayk Hall, one 


Fig. 4. Arjeri Cave Map — Russian Survey — app. 1985 


is actually standing near the top of the hall. 
The hall descends about 10 meters, and contains large, colorful columns, massive flowstone, and a myriad of 
formations. Overhead is voluminous white calcite coral — almost cloud-like in appearance. Continuing beyond 
Vayk Hall, there is a nearly 10 meter climbdown that can be negotiated with a handline. 
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Fig. 5. Arjeri Cave — Vrezh Nazaryan and Lara Fig. 6. Arjeri Cave — Lara Chavdarian in the Hall of 


Chavdarian in Vayk Hall - Republic of Armenia. 
Photo by C. Chavdarian 


Giants - Republic of Armenia. 
Photo by C. Chavdarian 


This then leads into a broad swath of massive columns, which we have appropriately labeled as the Hall of 


Giants. Continuing from this point leads one through a variety of colorful formations along the way, including 


nearly blood-red speleothems. Continuing on, the passage gradually slopes downward, and eventually leads 


to the “Lake” near the end of the cave’s known passage. On other trips into Arjeri Cave, we also explored 


some areas that were not present on the Russian map. There is much more exploration and detailed survey 


required of this cave, and it will surely extend the length of all known passages well beyond 4 kilometers. 


Karmir (Red) Cave: In 2007, from our base camp, we did a very steep hike of over 400 meters up the 


mountainside to Karmir Cave, located at an elevation of over 2.100 meters (Fig.7). The hike took over two 


hours and required careful negotiating of brush, loose talus, and some exposure along the way. Reaching the 


entrance requires care due to the final short exposed climb. 
Nearby is another cave — called Kiklop Cave — but the final 
30 meter climb to reach it has significant exposure and is 
also more challenging [4]. Time also prevented us from 
exploring this other cave. Karmir Cave has a large entrance 
chamber, allowing a group of cavers to congregate and also 
change in and out of caving attire. It is a conglomerate cave 
consisting of limestone and other mineralization, which 
provides the amazingly intense red color of much of the 
interior of the cave [5]. From the entrance chamber, there 
are several passages leading into the cave interior. The 
passages tend to intersect and loop back around. Well inside 
the cave we walked, crawled, climbed, and even traversed. 
The cave was wet and muddy, and permeated throughout 
with that remarkable red coloration. Although not highly 
decorated, there was some flowstone, stalactites, and red 
coralloids. We even observed some white (not red) 
moonmilk along the walls. There is also a 90 meter or so 
side passage in the cave, but neither our leader nor the rest 
of us could find it on this trip. The cave is well worth further 
exploration and is quite unique. 
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Fig. 7. Karmir (Red) Cave — Greg Chavdarian 
climbing in a red passage - Republic of Armenia. 
Photo by C. Chavdarian 


C. Chavdarian, S. Davtyan, S. Shahinyan 


Mageli Cave: Mageli Cave [4] is located along a gorge near the town of Areni, in Armenia’s wine country 
(Figs. 8,9,10). It sits at an elevation of nearly 1.100 meters and has over 2 km (1.3 miles) of passages. Mageli 
Cave is a classic example of a conglomerate cave, as it is a mixture of limestone and other minerals [6]. There 
is a map for the cave, but it is basically an overview, not a detailed representation of the cave. During our 
expeditions, this cave was explored on two different occasions — in 2007 and 2010. After entering the cave by 
crawling on hands and knees through the borehole entrance, one can resume walking through parts of the cave. 
Immediately, on entering the cave, the conglomerate nature of it is obvious, as the walls of the cave resemble 
coarse, pebbled or gravelly concrete. Some flowstone is observed, but, for the most part, the cave is devoid of 
decoration (and this is primarily due to the conglomerate composition of the cave). Not far inside the entrance, 
there exists a bat colony. Fortunately, because of the various passages, one can avoid passing near the bats. 
Exploring further into the cave, there are tall, narrow passages, boreholes, and a steep, slippery climbdown to 
a lower level (handline recommended). Well-inside the cave there is also a 3-meter long belly-crawl squeeze 
with an incredibly intense, cold wind blowing through it (similar to a venturi). We negotiated this squeeze and 
then ended up in a chamber that allowed us to stand. However, in examining the chamber, we were unable to 
locate any sizable opening or passage related to the heavy wind. After passing back out of the squeeze the way 
we came, and gradually working our way towards the cave entrance by a different route, we encountered long, 
booming passages and came upon an impressive, massive conglomerate block hanging down from above, one 
of the signature features of the cave. After exploring additional side passages, we exited the cave. In a 
subsequent trip to this cave in 2010, we entered through the main entrance and then proceeded to negotiate our 
way through passages, including boreholes, to the upper level of the cave (above the entrance). As we worked 
our way along this level, we encountered a 12 meter pit (but with no vertical gear, we did not drop the pit for 
further exploration). However, by passing through another nearby passage, the pit can be avoided, and one 
can, actually, exit the cave onto a narrow ledge above the main entrance (thus, an upper entrance). On this 
ledge, one has a spectacular view of the gorge below. A more detailed survey and subsequent map of this cave 
is certainly warranted. 


Mageli Cave 
M 1:1000 


Fig. 8. Mageli Cave Map — Survey — Smbat Davtyan 
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Fig. 9. Mageli Cave — Smbat Davtyan and Nyree Fig. 10. Mageli Cave — The massive conglomerate 


Abrahamian in borehole passages - Republic of Armenia. block with Seda Chavdarian - Republic of Armenia. 
Photo by C. Chavdarian Photo by S. Johnson 


Trchuneri Cave (Cave of the Birds): In 2007, we briefly explored Trchuneri Cave near the town of Areni. 
This cave is of particular note, as we observed evidence of an archeological excavation at the entrance of the 
cave [4]. Subsequently, in 2010 and 2011, significant archeological discoveries were reported at this dig, with 
artifacts dating back 6,000 years. Excavations continue to this day. 

Jerovank (Water Cave Church): From the base camp in 2007, we traveled lower in elevation to a gorge 
with a path, which lead to a very unique cave, known as Jerovank. After hiking through the gorge, we reached 
the small limestone cave, next to a running stream. The locals make a religious pilgrimage to this cave each 
year. A church altar was built inside the cave. Some of the nearby water actually permeates into the church, 
and the water can be seen along the back of the main chamber, and also in an alcove next to the church altar. 
Along one side of the church chamber a brick wall was built, further enclosing the cave and church from the 
outside. We observed flowstone along the natural walls of the cave. Due to centuries of turbulence in Armenia, 
church caves like this existed to provide protection for the congregation. 


Syunik Province, Republic of Armenia 


In 2011, we traveled to Syunik Province to photo-document the man-made caves of Khndzoresk and the 
man-made caves of Old Goris near the eastern border of the Armenia. 

Caves of Khndzoresk (Deep Gorge): The man-made caves of Khndzoresk are in an isolated and lightly 
populated mountainous rural region [7]. We hiked into this site. In centuries past, this was a thriving village, 
with man-made caves (homes) carved into the sandstone cliffs scattered all around (Fig.11). This even 
included a small church. It was inhabited up until the g 
1950s. The caves housed people, food supplies, and 
even livestock. All that now remains are the empty 
caves and the vacated church. Centuries ago, if one 
wanted a cave shelter or home, the local mason 
would fashion or carve a cave out of the soft rock. 
Over time, this resulted in a large, thriving cave 
village. The location of the caves rest on high 
ground, thus providing shelter in a strategic location 


— especially, important considering the enemies that 


Fig. 11. The man-made caves of Khndzoresk - Republic of 
Armenia. Photo by C. Chavdarian 


had invaded Armenia over the ages. 
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Caves of Old Goris: In the hills above the thriving 
Armenian town of Goris lie a series of rock spires 
(which look similar to “hoodoos”). The spires actually 
consist of many man-made caves, similar to those in 
Khndzoresk (Fig.12) [8]. This was the centuries-old 
village of Old Goris and it is extensive. The present 
town of Goris lies in a valley below Old Goris. One can 
access Old Goris, walk through it, and enter some of 
the carved sandstone caves. In some cases, where 
sandstone rock shelters may have already been present, 
the shelters were likely enlarged by local masons. The 


area is no longer inhabited. Once again it is important BATA aa 


man-made caves of Old Goris - Republic 
having provided protection to its former inhabitants. of Armenia. Photo by C. Chavdarian 


to note that Old Goris strategically sits on high ground, Fig. 12. The 


Tavush Province — Republic of Armenia 


In 2013, the focus was on the caves in the northeast section of Armenia, specifically in Tavush Province. 
We stayed in the town of Ijevan during our explorations in this province. 

Lastiver Cave (aka Anapat Cave): Lastiver Cave is perched near a cliff, in a mountainous region of 
Tavush Province, west of the town of Ijevan [9]. After driving up the mountainous terrain with our 4WD 
vehicle to a parking area, we then hiked about 4 km to the cave. The cave lies at an elevation just below 
1.200 meters. This natural cave is limestone-based. The cave has had some frequent visitation, and, as a 
result, is fairly devoid of decoration. The chambers in the cave range from about 10 to 100 meters deep to 
about 15 meters wide. In the past, one of the chambers was used as a crude church, dating back to the 12" 
century. In this chamber, and opposite the centuries-old altar, are a series of nearly two-meter tall, human 
wall carvings that appear to be just as old. However, it turns out that these carvings were actually created 
only about 80 years ago’, as explained by Smbat 
Davtyan. Why this was done is not clear. But it 
certainly has added to the mystique of the cave. In 
the largest room to the right of the church chamber, 
there is some actual flowstone remaining along the 
walls. However, there is breakage of the flowstone, 
and they are covered in dust and silt. This is now a 
fairly dry cave. In one of the other chambers we did 
observe two bats. However, there was no evidence 
that the cave houses a major bat roost. We explored 


various large and small chambers. Lastiver Cave is 


Fig. 13. Lastiver Cave — Smbat Davtyan and Vrezh 
Nazaryan in a large chamber with aged flowstone - 
(Fig. 13). Republic of Armenia. Photo by C. Chavdarian 


a natural cave that was also used as a church 


Large Grotto Cave, Pool Cave, and Crystal Cave: Hiking beyond Lastiver Cave, we gradually came 
upon a huge wall or cliff of karst at an elevation of about 1.200 meters, where we came upon a large cave 
shelter [10]. We named it “Large Grotto Cave”, as it consisted of a large entrance estimated to be 
approximately 45 meters in width. The cave extended back about 15 meters. There was breakdown at the 


> Speleological Center of Armenia 
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back of this large cave shelter. We looked for additional passage at the back of the cave but did not locate any. 
It is possible that there may be passage beyond the breakdown. 

We then hiked on along the karst wall and encountered another cave high up in the karst cliff. The way to 
this cave required some vertical rock climbing with exposure. We found a handline which had been placed 
there by others who had been to the cave. One member of our team — Vrezh Nazaryan — did the climb and went 
into the cave entrance. This cave measured about 7 meters in width at the entrance, and 10 meters deep. Further 
inside this cave (shelter), the width expanded to 35 meters. Inside the cave, someone had actually constructed 
a circular pool containing water. 

Hovk 5 and Hovk 1 Caves: The Hovk caves are high in elevation in a remote area of the mountains west 
of evan. There had been archeological excavations conducted 10 years earlier [4]. From our base in the 
town of Ijevan, we drove our 4 WD vehicle high up into that area. After parking off-road, we hiked up a very 
steep slope to the base of karst cliffs. We found Hovk 5 Cave at an elevation of approximately 2.380 meters. 
The cave is actually a rock shelter with dimensions of about 9 \ es : 
meters wide and 2 meters deep. We found evidence of an 
archeological dig. We then continued the hike near the karst 
cliff toward Hovk 1 Cave (Fig.14) [8]. During the hike, we 
noted limestone outcroppings scattered all along the way. Even 
more noteworthy was the lay of the land, which revealed an 
undulating landscape with sinks. As a result, there is great 
potential in this region for a potential cave system, or systems, 
waiting to be discovered. Once we reached the cave, we hiked 
up a slope and entered the elevated entrance, which sits at an 
elevation of over 2.050 meters. The first section of the entrance 
is about 0.5 meter wide and 3 meters long. This is followed by 
two large steps cut into the entrance — presumably by 
archeologists who excavated the cave entrance — which then 
lead into the main part of the cave which is approximately 3 
meters wide and over 21 meters in length. The main passage has 


a high ceiling which eventually pinches down to a marrow slot 


a 


Hovk I Cave — Vrezh Nazaryan 
clearly a very habitable cave, and strategically favorable as it in the high entrance - Republic of Armenia. 


sits up into the side of the karst cliff . Photo by C. Chavdarian 
Zrngan Cave (or Zerngan Cave): For our next 
destination, we traveled east of our base - the town of Ijevan 


at the back of the cave. There is little decoration. This was Fig. 14. 


- and up into the mountains towering above the town. The 
goal was to reach Zrngan Cave (Fig.15). at an elevation of 
approximately 1.850 meters [4]. The last part of our trip was 
off-road with our 4WD vehicle. This cave has a deep 
vertical entrance drop thought to be at least 45 meters in 
length. At the time of this trip, our colleague Smbat Davtyan 
told us that he knew of only one attempt at dropping down 
the entrance pit. This occurred about 30 years earlier. A 


group of non-cavers had lowered one of their people down 


Fig. 15. Zrngan Cave - C. Chavdarian on 
rope at the entrance drop. - Republic of Armenia. 
passage, but became frightened, and was hauled out of the Photo by V. Nazaryan 


to the bottom of the pit. This person actually explored some 
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cave. No one knows the true extent of the cave. Unfortunately, on the day we were there, C. Chavdarian was 
the only caver with extensive experience at single-rope technique. Although very tempting, the correct 
decision was made to not drop the pit — that is, not to do any solo caving for obvious safety reasons. Instead, 
he did a demonstration for the caving colleagues by rigging the entrance drop, gearing up, and only rappelling 
down a short distance, followed by a changeover and ascent out of the cave. An actual cave trip into Zrngan 
Cave may be planned for in a future expedition. 


Lori Province — Republic of Armenia 


After leaving Tavush Province, our final caving destination of the 2013 expedition was to a lava cave in 
the neighboring province of Lori, to the west. This is also a rugged mountainous area of Armenia that also has 
caves. 

Sanahin Lava Cave: Sanahin Lava Cave is the second largest lava cave (aka lava tube) in the Republic of 
Armenia (Fig.16), as noted by my colleague Smbat Davtyan [11]. The cave is located near the town of 
Alaverdi and is just off a main mountain road 
and on the side of a cliff, at an elevation of 
nearly 1.000 meters. The cave looks out over 
the town of Alaverdi and the Debed River in 
the spectacular and scenic gorge below [9]. 
This lava cave has about 80 meters of passage. 
There are three entrances to the cave, but only 
one is negotiable - and that is the main entrance 
nearest the road. To the left of the main 
entrance is a very small, squeeze entrance, 
which is not negotiable as a crawl, and farther 


to the left is a large, wide cliff entrance. The 


Fig. 16. Sanahin Lava Cave — Vrezh Nazaryan and Smbat 
Davtyan at the impressive cliff entrance, with a spectacular view 
vertical drop outside, and thus, very dangerous of the town of Alaverdi in the gorge below - Republic of 


if one is near this drop when standing inside Armenia. Photo by C. Chavdarian 


the cave and looking out this entrance. Over the years, people and livestock have inhabited the cave,but we 


cliff entrance has an approximate 30 meter 


saw no one there during our trip. The cave appeared to be abandoned. Inside, there is mostly walking passage. 
However, there is a wet passage of about 9 meters in length which required kneeling and crawling. We did 
observe some small secondary formations overhead in the cave inside this wet passage — small stalactites of 
about 2.5 to 5 cm in length, generated from solution deposition. Inside the main entrance and to the left is the 
main trunk passage. This passage curves and eventually ends at the wide cliff entrance with the sheer vertical 
drop. This cave is well worth exploring. 


Kotayk Province — Republic of Armenia 


In Kotayk Province sits the medieval monastery of Geghard. This church is noteworthy as it was originally 
created by boring by hand into a rocky mountainside. It is now a UNESCO World Heritage Site. The majority 
of the church sits inside the mountainside, thus, the interior walls and massive columns of the main chapel, 
built in 1215 A.D., are entirely carved out of rock. It is, in essence, Armenia’s most famous man-made cave 
church. Our team spent a day inside the church, marveling at it. 


The Nagorno-Karabakh Republic 


East of the Republic of Armenia lies the Nagorno-Karabagh Republic (NKR), home to an estimated 
150,000 Armenians. This present-day independent Armenian republic is actually part of historical Armenia 
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(and contains Armenian churches dating back centuries) but was tragically partitioned away and placed under 
the jurisdiction of Azerbaijan by the Soviet Union. Following the break-up of the Soviet Union in 1991, the 
Armenians of NKR fought for and regained their independence. This war with Azerbaijan ended in a cease- 
fire in 1994. Presently, one can travel through most areas of the Republic. Thus, in 2011, Vrezh Nazaryan, 
Smbat Davtyan, and Charles Chavdarian journeyed to NKR to ee a specific cave. 

The Caves of Tegh: On our way to 
NKR, we initially passed through the border 
village of Tegh (Fig.17), located just inside 
the Republic of Armenia, and east of the 
town of Goris. As we drove through, we 
looked off to our left, and there along a ridge, 
just below the main plain of the present-day 
village, we saw a number of large holes 
running along it. These were indeed caves, 
but it was not clear if they were all man- 
made. ome may have been natural caves that 
may have been enlarged to house families, 


and some may have been totally man-made. { i 
We did not stop there, but we did capture Fig. 17. The Caves of Tegh - Caves below the village - Republic 
some photos. of Armenia. Photo by C. Chavdarian 

Azokh Cave: Azokh Cave, located in the southern region of NKR, was the goal of our trip. Lying on a 
hillside in NKR’s Hadrut province, and overlooking the village of Azokh, is Azokh Cave (Fig.18). This cave 
is of archeological significance, as a Pleistocene, pre-Neanderthal mandible fragment was discovered there in 
1968 [12]. This led to a series of subsequent excavations up to the present. A number of artifacts have been 
found, but the excavations are essentially all at the three entrances (not deep inside the cave), with most of the 
excavations at the large main entrance (which is a high and wide vertical slot entrance). It is known that the 
cave consists of roughly three chambers, al 
and three cave entrances, and has roughly 
180 meters of cave passages. In 
researching this cave prior to the 
expedition, we found very little 
information regarding the actual interior of 
the cave (beyond the main entrance), and 
hardly any photographs. There is only a 
very rudimentary hand-drawn outline map 
of the cave [5]. It appeared that very few 
individuals ventured beyond the main 


entrance and deep into the cave. Once we 
passed through the main entrance chamber Fig. 18. Azokh Cave - Smbat Davtyan with a column and 
and into the next chamber, we discovered decoration - Nagorno-Karabakh Republic. Photo by C. Chavdarian 

why there was so little information. In the second chamber and beyond, one is subjected to massive amounts 
of flying bats, flies,and guano (and its odor), which can be overwhelming. This was a shock, but it also 
answered the riddle of this cave, as to why so little was known of it. A non-caver would not have lasted longer 
than 10 minutes inside the cave. However, we cavers were determined to explore and photo-document 
the cave. We managed to remain in the cave for about 1.5 hours. We gingerly moved through the narrow 


passages connecting the chambers, either upright or by crouching, and were careful not to fall into the slippery, 
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deep guano. The cave is mostly devoid of decoration however, we were pleased to find two large limestone 
columns — one of which was rather impressive with surrounding flowstone. We also discovered an incredibly 
large mound of guano taller than us — reminiscent of the guano mound in the Planet Earth television series. As 
tough as this trip was, all caves deserve to be explored, regardless of the conditions one may encounter. We 
accomplished the objective, and now have a clear understanding of the nature of this cave. 


Discussion 


Inside one of the caves, someone had actually constructed a circular pool containing water. Stones lined 
the pool. As a result, we named the cave “Pool Cave”. It is not clear the purpose of the pool. It may have 
provided a water source for the more contemporary inhabitants, or possibly a source for bathing. 

We then continued our hike along the bottom of the karst wall and soon came to another cave at an elevation 
of approximately 1.200 meters. The walk-in entrance was approximately 2.5 meters high and 2 meters wide. 
The cave extended back over 20 meters, and we could walk upright through the entire length of it. There was 
no noticeable decoration, and a fair amount of breakdown. It is possible that there had been an excavation 
inside this cave. The notable feature was that there was a pile of crystalline calcite near the entrance, which 
had been partially scavenged by others. Because of this, we simply named the cave 
“Crystal Cave”. 

The medieval monastery of Geghard in Kotayk province is, in essence, Armenia’s most famous man-made 
cave church. The original name of the monastery was Ayrivank, which means “The Cave Monastery” or “The 
Monastery of the Cave”. It was later called Geghard, which refers to the spear which wounded Christ during 
the crucifixion, and was subsequently alleged to have been brought to this church for storage (the spear 
currently resides in the Holy See of Echmiadzin, near the capital of Yerevan, Armenia). In any event, Geghard 
is a must-see for any visitor to Armenia, including cavers. 


Conclusion 


Over a period of six years, with four expeditions to Armenia in 2007, 2010, 2011, and 2013, we explored 
and photo-documented a number of natural and man-made caves — limestone, conglomerate, sandstone, and 
lava. We experienced the color and beauty of Armenia’s natural caves, the intriguing man-made cave villages, 
and the reverent use of caves as chapels and churches. Caves have been an integral part of the landscape and 
culture of Armenia and the Armenians. Through these expeditions, one cannot help but gain a greater and 
enduring respect for this ancient land and its people. 
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future existence also became reliant on the sufficient availability of energy ~ 

and its sources. Nevertheless, energy consumption is becoming a major 

concern due to the increase in usage rates and the lack of sufficient renewable Received: 16.02.2023 

sources of energy. Nowadays, the technological advancement that lies in the ; 

use of photovoltaic panels (PV panels) can help in generating energy and Revised: 10.03.2023 
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multiple factors that can ensure the highest potential generation of energy. 
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the panel. Modifying the titling angle can play a significant role in generating 
high rates of power depending on the measurements and calculations. The 
aim of this paper was to investigate the efficient titling angle that can be used 
for PV panels that are installed in Al-Sherouk City in Egypt. The methodology 
involved the implementation of an experimental investigation that is based on 
position 2 similar PV panels to estimate the power generation over a period 
of 2 days from 9 am to 4 pm. The findings have shown that the theoretical and 
experimental results were similar, and the optimum tilting angle was 
determined to range between 54.7 degrees and 8.16 degrees. The study 
demonstrated the differences between an adjusted and a fixed angle, and the 
variable titled angle can generate more energy than the fixed one. This paper 
contributes to the body of knowledge by presenting the significance of a 
variable tilting angle to generate more power than relying on a fixed angle as 
demonstrated by many previous studies. 


This work is licensed under a Creative 
Commons Attribution-NonCommercial 
4.0 International License 


Keywords: solar energy, photovoltaic panels, renewable energy, solar 
energy in Egypt, optimum tilt angle 


Introduction 


People should pay more attention to renewable energy due to the lack of sufficient energy sources and the 
positive impact of using a renewable source of energy. It is important for developing and developed countries 
to consider finding new solutions to increase in prices of oil, increase in electricity demands, high rate of 
greenhouse gas emissions, and even global warming [4]. Solar energy is categories as a renewable source of 
energy which has been developing to an extent scale due to limitations in energy transmission [2]. Normally, 
there are multiple benefits of solar energy over the use of fossil fuels such as reduced carbon emissions, cleaner 
air, and can generate power over a long period of time. Hence, due to the huge increase in electricity 
consumption, researchers are more concerned about developing solar energy developments with excellent 
efficiency, less environmental pollution, and with proper investments cost [2]. 

The overall global energy consumption in 2014 was recorded to be around 160310 million MWh, and it 
was estimated that this value is expected to keep on increasing to reach around 240318 million MWh by 2040 
[6]. It was recorded in 2010 that the electrical generation systems that are based on renewable energy increased 
by 20% from the overall electricity generation, and this value will keep on increasing until it reaches 31% by 
2035 [3]. The use of renewable energy can help in generating around 57% of global electricity supply by 2025 
[3]. Eventually, the dependency on solar energy for electricity production is growing around the globe. Hence, 
solar energy production should be developed and promoted due to the negative impact of conventional energy 
production [1]. During the past few years, a huge amount of investment was focused on the improvement in 
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solar energy production which helped countries to reach better technological advancement and cost-effective 
production of energy [1]. 

Shu et al [5] investigated the optimum tilt angle for solar panels in Kitakyushu City. There are multiple 
variables that were considered in the experiment including the sensitivity of optimum tilting angle, radiation 
rates, reflection rates, and declination of solar. The results demonstrated that the 35 degrees tilting angle had 
the most intensity of radiation between December to November of the following year. Another important factor 
to be considered is that if the latitude gets higher than the optimal angle might increase. 

Tlijani et al [7] conducted a study regarding the optimization of tilt angle of solar panels that are located in 
Tunisia. The research focused on the comparison between the experimental and theoretical results using Matlab 
simulation analysis. The angles used in the study were 0, 30, 45, 60, and 90 degrees and the panel was 
positioned once in the west, then changed to south and east. The first panel had dimensions of 370*295*15 
mm and was placed at a temperature of 25 degrees Celsius. The second one had a dimension of 215*225*18 
mm and was located in the same way as the first one. The findings of this study demonstrated that climate 
conditions and the location of the solar panel might impact the total power generated by the system. This paper 
discussed the potential optimum tilt angle for PV to generate optimum power and energy. 


Materials and Methods 
Procedure of Calculations (Table 1) 


1. The first step taken in the calculations was determining the deflection angle as shown in the equation 
below: 


() 


284 + *) 
365 


Deflection angle equation 6 = 23.45 sin (3 60°: 


2. Bis then calculated using the following equation 


360 


pega 
OY Vane 


3. True solar time is then measured using this equation 


E = 229.2 - 0.000075 + 0.001868 - cos B — 0.032077 - sinB — 0.014615 - cos2B — 0.04089 - sin 2B) 


4. After calculating the B, the solar time (ST) is then calculated, ST is the time where the measurements 
are taken (14:00pm). This time was basically chosen because the radiations from the sun were 
considered the highest during the entire day time. Is, is the standard longitude of time zone, and I loc is 
the longitude of the location. Io: of the chosen location “Al-Sherouk City” is taken as 31.610584. 


Solartime — ST = 4° (Is¢ — ioc) + E 
5. Hour angle is then measured using the following equation, where the t, is known as the solar time. 
Hour angle w = 15(t, — 12) (2) 


6. Zenith angle is then measured as follows, where ® is the latitude and taken as 30.1187 according to Al- 
Sherouk City location. 


Zero angle cos 6, = cosd:cos@:cosw + sind « sing (3) 
7. The solar elevation angle is then measured using this equation. 


Calculating the zero elevation angle a=90-86, (4) 
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8. Solar Azimuth angle is the next to be calculated. 

Solar azimuth angle Y; = cos~*[(sin(a) « sin(@) — sin(5)/cos(a) - cos(@) (5) 
9. Incidence angle is then measured using the following equation. 

Incidence angle 6 = cos~*[cos(6) « cos(@) « cos(w) + sin(6) « sin(@) (6) 
10. Then the tilt angle is measured as follows. 


Calculation for the tilt angle B=|G-6| (7) 


Table 1. Measured calculation of all above parameters 


Days / 6 B E Solar W Cos | Zenith Solar Without Solar B 
parameters time zenith | angle elevation (cos) azimuth 
angle angle (a) inverse | angle 0's) 
1 -23.01| 0.00 |-2.90 | 14.16 | 32.34 | 0.46 62.37 27.63 0.61 52.45 54.62 
2 -22.93| 0.99 |-3.35 | 14.16 | 32.45 | 0.46 62.35 27.65 0.61 52.53 54.54 
3 -22.84| 1.97 |-3.39 | 14.17 | 32.56 | 0.46 62.33 27.67 0.61 52.60 54.45 


Calculating the solar irradiance 


The theoretical calculations are based on previously obtained results of using solar panels in the same 
locations investigated in the experimental part of this study. The below are the data collected from specific 
dates in 2017 (Table 2). 


Table 2. Calculation of Solar Irradiance of tilted angle "S" 


Date A B C Zenith Tilt Gbn Gd G rb F > Gt 
angle | angle Zs 
23 
Os 
2 
17/1/2017 | 1232.46 | 0.1416 | 0.057 61.31 52.23 | 917.46 | 52.93 | 493.38 1.02 0.988 0.2 523.4 
16/2/2017 | 1216.58 | 0.1436 | 0.059 55.19 | 44.57 | 945.87 | 56.44 | 596.44 1.01 0.991 0.2 622.8 
16/3/2017 | 1190.18 | 0.1538 | 0.009 45.99 34.03 | 953.76 | 52.04 | 714.69 1.00 0.994 0.2 733.2 


15/4/2017 | 1144.68 | 0.1753 | 0.091 35.69 22.2 | 922.41 | 84.83 | 834.1 0.99 0.989 0.2 844.60 


15/5/2017 1109.4 | 0.1928 | 0.116 28.87 | 12.82 | 890.17 | 103.43 | 882.97 | 0.98 0.986 0.2 884.13 


11/6/2017 | 1092.84 | 0.2020 | 0.129 27.28 8.52 | 870.56 113 886.74 | 0.98 0.983 0.2 884.09 


17/7/2017 1085.4 | 0.2067 | 0.135 29.76 | 10.43 | 855.39 | 116.10 | 858.68 | 0.98 0.981 0.2 859.78 


16/8/2017 | 1103.45 | 0.1967 | 0.124 33.9 18.16 | 865.15 | 107.50 | 825.59 | 0.99 0.983 0.2 834.43 


15/9/2017 | 1142.48 | 0.1816 | 0.097 40.06 | 29.39 | 901.11 | 88.13 | 777.82 | 0.99 0.987 0.2 795.64 


15/10/2017 | 1183.8 | 0.1634 | 0.076 48.3 41.21 | 925.98 | 71.11 | 687.10 1.00 0.989 0.2 713.75 


14/11/2017 | 1213.68 | 0.1514 | 0.052 56.51 | 50.52 | 922.30 | 60.18 | 569.10 1.01 0.988 0.2 601.07 


10/12/2017 | 1228.23 | 0.1445 | 0.059 61.24 | 54.66 | 909.47 | 53.84 | 491.42 1.02 0.987 0.2 524.11 
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Table 3 represents the average day for each month in order to measure the solar irradiance in this 


specific day. 
Table 3. Average day for each month 
penis Wiorih For average day of month 

day of month Date n 6 
January i 17 17 -20.9 
February 31 +1 16 47 -13 
March 59 +i 16 75 -2.4 
April 90 +i 15 105 9.4 
May 120 +i 15 135 18.8 
June 151 +i 11 162 23.1 
July 181 +i 17 198 21.2 
August 212 +i 16 228 13.5 
September 243 +1 15 258 2.2 
October 273 +1 15 288 -9.6 
November 304 +1 14 318 -18.9 
December 33441 10 344 -23 


Calculating Power Output and Efficiency of Panel: 


The total power output had to be measured for the panel during the experiment for the whole month 
according to the data obtained from the British University in Egypt as shown in Table 4. 


Table 4. Measurements of power for each month 


Jan | Feb | Mar | Apr | May Jun Jul Aug Sep Oct Nov | Dec 
os Optimum | 57.5 | 53.1 | 584] 55 | 564 | 546 | 560 | 55.5 | 53.9 | 55.9 | 53.7 | 55.6 
(KWH/month) 
Power : 
Tilt 30 52.4 | 50.9 | 58.2 | 54.3 | 53.2 | 50.0 | 52.1 | 54.0 | 53.9 | 546 | 49.6 | 49.6 
(KWH/month) 


Maximum efficiency equation 


Pnax (maximum power output) 


ga OBERT SECA) = (E (incident radition flux) - A.(area of collector) (8) 


Experimental Procedure 


The experimental procedure that was followed in this 
research focused on comparing the generated results between 
the theoretical analysis and experimental results and indicate 
any variations in the generated results (Fig.1). The experiment 
included the use of 2 identical panels that are directed to the 
south direction using multiple tilt angels that ranged from 0 


degrees, 15, 30, and until 45 degrees. These panels were | a 
positioned on the top of the British University in Egypt. Fig. 1. Used PV panels in the experiment 
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The main objective was to indicate the solar radiation and generated power from the panels during an annual 


measurement where the size of the PV module used was 1.65 - 0.953 meters and the overall area is 1.57 m2. In 


addition, the maximum power within the panel was 280W. The type of the used PV panel is Monocrystalline 


Silicon Module. The field observations were decided to be undertaken during 2 various days which are 
24/5/2018, and 26/5/2018. The PV panel was installed and kept on recording starting from 9 am until 4 pm, 
while the readings are taken each half an hour. The field procedure consisted of the following steps: 


1. The first step was to properly clean the panels in order to remove any sand particles, dust, or any objects 


that might affect the collection of solar radiation by the panel. 


2. Then, the panels were adjusted into 2 different angles at which to determine the optimum angle that 


could generate the highest solar radiation. 


3. PV system analyzer was adjusted to include the data and specifications of the panel used. 


Experimental Results 


The results demonstrated that 
15 degrees tilt angle generated 
the highest potential power 
output from the panel as shown 
in Fig.2. 

Fig. 2 shows the experimental 
readings that were taken during 
the 2 days of field observations 
starting from 9am until 4pm. 
Certain radiations were chosen 
after the experiment in order to 
estimate the power energy while 
the most optimum angle will 
usually yield the highest power. 
In the theoretical results the 
optimum angle was estimated to 
be from 10 to 11 degrees as 
shown in Fig. 3. The highest 
angles that provided most power 
were recorded in 0, 15, 30, 45 
degrees respectively. Therefore, 
the experimental and theoretical 
results were quite close 
according to the generated 
optimum tilt angle (Fig.3). 


400 


experimental reading 


450 500 550 600 650 700 750 


Radiation(W/m2} 


——angleO —#—angle 15 angle 30 —®—angle 45 


Fig. 2. Experimental readings for the different tilt angles and power outputs 


optimum tilt angle 


Days 


Fig. 3. Tilt angle for May (theoretical observations) 


Difference between Theoretical and Experimental Results 


Radiation Results 


35 


The rate of radiation was measured for all angles and compared with the theoretical values as shown in 


Figs. 4,5,6 and 7. 
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Fig. 4. Difference between experimental and theoretical radiation at 0 degrees tilt angle 


The theoretical radiation was optimum during the mid-day and kept on reducing by time. The experimental 
line was also similar but with less solar radiation during most of the day except during mid-day. 


angle 15 degree angle 45 degree 


@ theoretical radiation — experimental radiation 


r 
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radiation 


time time 


Fig. 5. Difference between experimental and theoretical __ Fig. 6. difference between experimental and theoretical 
radiation at 15 degrees tilt angle radiation at 45 degrees tilt angle 


If the tilting angle was rotated to 15 degrees, the theoretical readings were quite similar to the experimental 
ones as shown in Fig. 5, and the 30 degrees tilt angle also yielded similar results to this one. 

The 45 degrees tilting angle scored different values than the theoretical numbers as indicated in Fig. 7. 
According to the previous figures, the results were quite close and identical in terms of solar radiation which 
proves the credibility and validity of the generated results. However, it was observed that most of the 
experimental readings suffered from potential weather conditions such as dust. The highest solar radiation was 
estimated to occur at the solar time which is around 12 pm mid-day. Another measurement that was considered 
in the experiment was the power generated by the panel. The highest variation that was observed between the 
theoretical and experimental results was estimated to occur at 11 am. The differences between the generated 
power between the experimental and theoretical results was around 73 Watt. The following are the results of 
each tilt angle and the generated power: 

e In 30 degrees angle the power variation was 66 Watt at 11 Am. 
e In 15 degrees angle the power variation was 73 Watt at 12 Pm. 
e In 45 degrees angle the power variation was 41 Watt at 10 Am. 
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The results have shown that the major 


Energy output 


variations between theoretical and 
experimental results were observed between 


10 am and 12 pm due to the efficiency of the 3 2 
panel. It must be taken into account that the 2” | | i | | | i | 
theoretical observations are taken with the a 
optimum panel efficiency, while the pe 
jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 


experimental ones might vary according to 


Month 


the properties and conditions of the solar moptimum mtitt30 
panel. Fig. 7. Power vs months output 

Fig. 7 demonstrates the generated power by solar panels for each month. It helps in understanding the 
generated power from the optimum tilt angle during the entire year. The blue bar shows the optimum energy 
output which keeps on changing during the year. The orange bars indicate the energy output for the solar panel 
if the tilt angle is set at 30 degrees. The used tilt angles for these solar panels were 
50,45,34,23,13,9,11,18,30,41,50,55 respectively. The results indicated that this program can be used to 
determine the critical tilt angle, solar radiation, and most possible power through changing the latitude and 
longitude of certain points any location. 


Parameters that Affect Theoretical and Experimental Results 


There are multiple factors at which play a major role in introducing variations between the theoretical and 

experimental results such as the following: 

1. Weather conditions due to the presence of clouds that covers the sun most of the day. 

2. The radiation sensor is not accurately fixed in the system at which can influence the final readings. 

3. The ground reflectivity factor is not quite accurate at which can generate the highest power possible. 

4. Calculating the A, B, and C factors is not quite effective as it might contain any errors especially during 
any calculations. 

5. Assuming that the rb ratio is zero when the incidence angle is higher than 90 degrees is not quite effective 
and can generate multiple errors. 

6. The efficiency depends mainly on the solar radiation and generated power which ranged from 7.9% until 
12.85% in the experimental results, while in the theoretical calculations the efficiency kept on the same 
value of 17.8% which is the highest efficiency possible. 


Conclusion 


Renewable energy is the source of electricity generation in the future and consistent research and 
development is required to enhance the current knowledge and understanding. One of the most common 
sources of renewable energy is solar panels and further studies are essential to better understand the design and 
efficiency of this system. Future researchers that are looking to apply the same conducted theory must take 
into account the results obtained in this research. It is important to take solar radiation results and 
measurements over a long run as it can be affected by the weather conditions or any other external factors that 
would prevent the generation of accurate results. Designing the solar panel to have a fixed inclined angle might 
help in generating more energy, and studies also recommended developing panels that can flexibly rotate 
during the day to generate optimum energies. The study investigated the theoretical and experimental results 
of power generation of solar panels in Al-Sherouk City in Egypt. The findings have shown that the theoretical 
and experimental results were similar, and the optimum tilting angle was determined to range between 54.7 
degrees and 8.16 degrees. The study helped in providing the significance of the tilt angle for power and energy 
generation of solar panels. Future studies are expected to investigate similar properties in other regions and in 
different climate conditions. 
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Recommendations 


There are various factors to be considered by future studies in order to develop on the findings that were 


reached in this study including the below: 


1. 


Changing the tilt angle can generate the highest possible solar radiation and power, thus it must be 
properly applied and taken into account, but at the same time the location of the panel and weather 
conditions might play a major role in finding the optimum angle. 

Ensure the consistent change of tilt angle monthly to generate the highest energy possible. These angles 
are 50,45 ,34,23,13,9,11,18,30,41,50,55 respectively in El-Sherouk City. 
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Abstract: Artificial cavities, i.e. man-made structures excavated within Bixio Roberto* 

rock masses in the mountains, below the ground, or in the subsoil of urban E-mail: roberto_bixio @ yahoo. it 

areas, are typically distinguished based upon the epoch of realization and 

the function for which they were originally used. They can be ranked into 

the following types, in turn divided in sub-classes: hydraulic works, 

dwelling works, worship works, war works, mining works, transit way Received: 14.03.2023 

works and others. The above criteria are essential for establishing a Revised: 04.04.2023 

general common line aimed at providing optimal elements for , ae 

cataloguing and comparing subterranean features, which may favor the Accepted: 20.04.2023 

creation of databases functional to knowledge, protection and 

enhancement of the hypogean works. In addition, there is another useful 

aspect for studying the origin and evolution of underground structures 

that takes into account their implementation modalities. The National 

Commission on Artificial Cavities of the Italian Speleological Society has © The Author(s) 2022 

identified, according to its experience in the field and in function of the 

construction techniques, six general categories of underground works: 

cavities dug in the subsoil, cavities built in the subsoil, cavities obtained 

by re-cover, anomalous artificial cavities, mixed artificial cavities and This work is licensed under a Creative 
. . : pus Commons Attribution-NonCommercial 

natural caves modified by men (anthropized caves). In this contribution AW inierational License 

we will discuss the specific details of each category, thus extending the 

concept of rupestrian heritage, usually confined to temples or dwellings 

carved in the rock, to a culture of building in "negative" that finds larger 

and more diversified evidences. 

Keywords: Rupestrian Works, Artificial Cavities, Categories, Typologies, 

Underground, Construction Techniques. 


Introduction 

And nowadays, man-made structures excavated within rock masses in the mountains, below the ground, or 
in the subsoil of urban areas (artificial grottos, cavities), do not cease to attract the researchers’ attention, 
although they are sufficiently studied and classified based on the epoch of realization and the function for 
which they were originally used. 

For those who study the origin and evolution of underground structures, accurate records of the ways of 
their implementation become very important in order to correctly understand the scientific problem, its 
significance and direct the study to certain results. As such, six general categories of underground work defined 
by the National Commission on Artificial Cavities of the Italian Speleological Society, according to its 
experience in the field and in function of the construction techniques are considered, which, in their important 
details, expand the hitherto existing conception of the rupestrian heritage. 


Materials and Methods 


During the preliminary research, data in the scientific literature was examined. Materials about the cavities 
located in Armenia were developed from publications in international conferences. The necessary sections of 
geological, stratigraphic, geomorphological, hydrogeological, and geophysical maps were collected and 
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developed. During the fieldwork of the scientific expeditions, a comparative analysis of these maps was carried 
out using data gathered from each site. Then the main works began, such as identifying the artificial caves, 
determining their location in the rock mass, comparing the conformation of the cavity elements with 
appropriate topographical surveys, and drawing the relevant conclusions. 


From “rupestrian works” to “artificial cavities” 


Usually, when the “rupestrian works” topic is dealt with, one thinks immediately to structures dug by men 
in the rock-faces in order to obtain underground spaces to inhabit, with related sites for productive activities 
(warehouses, stables, mills, dovecotes, etc.), or to spaces used for worship, both for liturgical purposes 
(temples, churches, etc.) and for the burial (tombs of various types), according to a functional types subdivision 
[1-4]. 

We find striking examples in various parts of the world, as the rupestrian settlements in southern Italy, the 
most famous of which are the Sassi of Matera, consisting of thousands of dwelling units. However, not less 
important are the many sites excavated in the gorges (locally named gravine [5]) of the Apulian-Lucan area, 
counting almost 600 churches, nearly the same number estimated for the hundreds of rupestrian settlements of 
Cappadocia, in central Turkey, and others in Armenia. Equally well known are the Dogon cliff villages of 
Mali, the Buddhist “caves” of China and India, the rupestrian city of Petra, in Jordan, the Pueblo villages in 
Colorado, and many others in different parts of the world. 

The use of exploration techniques derived from the experiences of cave progression has allowed, since the 
1960s, the finding of underground man-made works, less visible than those mentioned above, more dangerous 
to explore, and more difficult to be documented. This approach has produced a quantitative increase in 
knowledge of underground structures and, especially, has extended the investigations to issues not much 
considered before about the use by man of underground structures. This union, which began in prehistoric 
times with occupation of the caves by man, later evolved over the millennia with surprising and ingenious 
works. These are not limited to the already mentioned cavities intended for dwelling or worship, but include 
structures for transit, hydraulic engineering, mining, war, and related sub-classes (see Table at the end of the 
article). 

Furthermore, the research field has expanded from the works excavated by man in the rocky outcrops or 
below the countryside level to those in the subsoil of urban areas. 

It was therefore necessary to go beyond the idea of “rupestrian work’, as defined above, replacing it with 
the most extensive “artificial cavity,” or “anthropic cavity”, in its turn complementary to that of a natural 
cavity, or cave. We can define artificial cavity as a space created by man in the subsoil, in broad sense, which 
implies an idea of a "negative" construction culture, as an alternative to the outer buildings at the surface, but 
also to the hypogean environment produced by meteoric agents and geological phenomena. 

This approach completely defines the relationship between man and the underground world, taking into 
account not only the many variable purposes - identified in the types above described — but also the different 
construction ways according to morphological, lithological and urban characteristics of the environment in 
which the structures have been made over time, also including those structures which were not excavated, but 
share many similarities with hypogean places. 


Categories of construction techniques 


The Commissione Nazionale Cavita Artificiali (CNCA), i.e. National Commission on Artificial Cavities of 
the Societa Speleologica Italiana (Italian Speleological Society), according to the exploration experience of 
its researchers, has developed a catalogue of artificial cavities based upon the construction techniques, 
identifying six categories, according to their intended purpose [6]. In the Register of Artificial Cavities, 
compiled by the Commission [2,7-9], not only the strictly rupestrian works are therefore included, but all those 
structures built or dug by man in the subsoil and, sometime, in the above ground, according to the following 
criteria. 
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Cavities dug in the subsoil 


These are obtained exclusively by removal of stone materials (rocks), and can be divided into two 
groups [10]. 

a) Rupestrian (or rock-cut) structures, strictly 
speaking. They consist of spaces (rooms, tunnels, 
shafts) dug by man above the ground surface, in the 
outermost portion of rocky towers, pinnacles, cliffs, 
canyons, slopes. These are also defined “cliff 
cavities” [11]. 

Generally they can have very long horizontal 
development, even kilometres, on a single level 
(linear cavities/settlements), or on a series of 
stepped levels (terraced cavities/settlements), or on 
levels superimposed over the same vertical wall 


(wall cavities/settlements). When dug _ inside : oe ties Mo: =e 
individual pinnacles they are called "cone cavities" Fig. 1. Cliff rock-cut village of Hasankeyf on the Tigris 
(Fig. 1,2). River, in southeasternTurkey (photo M. Traverso) 


b) Underground structures, dug in depth (deep layer), under the ground level (mesas or plateau areas) or in 
the inner part of rock mountains (butte and other ridges). In this case, too, we can have networks extending on 
a single horizontal level (Figs. 3,4), or networks descending in the subsoil for tens of meters on superimposed 
levels (Figs. 5,6). 


CLIFF ROCK-CUT VILLAGES - VILLAGGI RUPESTRI IN FALESIA 


linear settlement terraced settlement __ wall settlement _ cone dwellings _ 
insediamento lineare insediamento a gradoni insediamento a parete insediamento a coni 
= F 
<r pit graves 
~ 
sz Yaa) Pigeon-loft 
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ashi () Br ehaiee dwellings 
or stores 
De 7 trench - stable 0" I 


; 0 
an SE, apiary on 


(winery) (? welling 
coy ; eh OB bese hapel 
oben, ad ae horizontal shelter oil mill 2 nt 
silo ; crusher orchards 
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ire a tunnel 
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Fig. 2. Exemplification of different cliff rock-cut village models in Cappadocia, 
central Turkey (drawing R. Bixio) 
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Fig. 3. Scheme of the "horizontal" underground shelter Fig. 4. Artificial tunnel under the site of Ani, 
of Filiktepe in Cappadocia, central Turkey ancient capital of Armenian kingdom, now in - 


(drawing R. Bixio) eastern Turkey (photo R. Bixio) 


“ae 


Fig. 5. Scheme of the "vertical" underground settlement Fig. 6. Underground settlement of Derinkuyu in 
of Derinkuyu, in Cappadocia, central Turkey Cappadocia. Well at depth of about 40 m 
(drawing R. Bixio) (photo G. Bologna) 


Note: the distinction between rock-cut cavities and underground cavities is not always clear, so that often 
the terms can be used interchangeably. 


Cavities constructed in the subsoil PHASE 1 PHASE 2 


The underground spaces belonging to this category are 
those obtained with masonry works created to define terrace 
volumes produced as a result of excavation of the subsoil, Fate ot ekavatcn 
through two techniques. 


c) Tunnel excavation technique (tunnelling). Removal 


of the rock is carried out entirely underground. The rooms 
are then coated with different masonry techniques. Coating 


can interest only part of the excavation. 


PHASE 3 


d) Trench excavation technique. It is realized with a 
open air excavation, followed by total or partial coating of 
the walls, building of the vault, and finally re-covering 
(Fig. 7). It is a technique very useful at not great depths: 
generally, it is simpler, faster, and cheaper than the tunnel denouing 
excavation technique. It can also be used in clayey soils. 


In some cases, the walls are not coated, and the only 


built part is the cover, which can be with flat ceiling, 
obtained by laying stone slabs or concrete slabs, or barrel- Fig .7. Reconstruction of the excavation technique 


like ceiling, obtained using various materials such as of the trench of Ahlat, lake Van, eastern Turkey 
ashlars, bricks, or concrete. (drawing R. Bixio) 
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The underground aqueduct of Gravina in Puglia 
(Southern Italy), for example, has been realized with 
both the tunnelling technique, without any coating, 
and the trench excavation technique, with partial 
coating and barrel covering by tuff blocks (Fig. 8) 
[12,13]. 


Mixed artificial cavities 


These are works dug to reach, extend or modify 
substantially a natural cave. The artificial part can be 
carried out with one of the methods described above. 

We found a significant example of mixed cavity 
in the archaeological site of Troy (Turkey). This is a 
work for the low city’s water supply, attributed to 
the 3 millennium BC and used until the Byzantine 
period [14]. The initial part, accessible from the 
outside, is a natural tunnel. The artificial 
part consists of a long tunnel, almost straight, led 


solid rock 


from an internal point of the cave to the base of some 
ascending shafts/wells, dug into the body of the rock. 
The excavation of the tunnel was made using the 
opposite fronts technique, that is, by two teams 


digging one towards the other, most likely 


simultaneously: this is evidenced at the junction Figs 8 Undeveromnd anueduerar Gravnain Piha 


point, typically identified because of the change in (soythern Italy). Mixed excavation technique in solid 
direction and the related blind appendix due to a rock: tunnel and trench, with partial coating and 
slight alignment error (Fig. 9). cover in stone ashlars (photo M. Traverso) 
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Fig. 9. Plan of Troia water supplying tunnel, mixed cavity (R. Bixio, elaboration after M. Korfman 2003) 
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Anthropized caves 


We define “anthropized caves” the natural caves 
that have undergone limited human interventions. 
They represent the boundary line between the natural 
cavities, produced by weathering and geological 
phenomena, and the artificial or anthropic cavities, 
entirely man-made in the subsoil. 

In general, these are wide but not very extensive 
caverns, in which man has built masonry structures, 
for dwelling and/or for worship, sometimes 
supplemented by small digging actions. 

The best known examples are, for the first type, the 
“Pueblo” villages built by American Indians in 
Arizona, Colorado and New Mexico between the 12th 
and 14th centuries (Fig. 10). cae 

The second type is represented by the sanctuary- Fig. 10. Antréphized 


ke 


cave: Pueblo at Mesa Verde in 
caves. We have examples in Italy (Santa Lucia, Colorado, USA (photo G. Stalteri) 


Toirano), France (La Saint Baume, near Var), Turkey (Sumela, near Trabzon), to cite a few, and in many other 
places in the Mediterranean basin, where there are churches and monasteries from the early centuries of 
Christianity, built inside karst caves. 


Non-excavated artificial cavities 


Finally, we describe the “re-covered cavities” and the “anomalous cavities” that are quite different from 
the rupestrian works, as previously defined, because they do not contemplate digging works to obtain spaces 
within the rock mass. However, they are included in the classification of the Commissione Nazionale Cavita 
Artificiali and inserted in the related register, as there is no doubt that, by their nature, they fall within the 
categories of anthropic cavities. 


Cavities obtained by re-covering 


Often the human activity on surface, particularly in BA 2 

~~ “concrete slabs ~~. -—— 
urban areas, has produced the overlap, the burial and : 
the embedding of natural or artificial spaces originally 
not located in underground spaces. 

For example, the fifty-two streams crossing the city 
of Genoa (Italy), in medieval times flowed in sub- 
aerial beds [15,16]. The need, with urban growth, to 
obtain new spaces for the city, has caused over the 
centuries their progressive coverage (Fig. 11), almost 


or Si ines “ 
always coincident with road axis, producing the Pemuralriverbed CS. A 


«. 7, , “ P 


incorporation of existing structures (bridges, dikes, Fig, 11. Rio Groppo, Liguria, Italy. Tunnel obtained 
masonry banks, remains of buildings, etc.). with the cover of the river bed (photo M. Traverso) 

For the sake of completeness, we would like to remind that there are also natural underground waterways 

(karst rivers) throughout the world. In this case, they are considered as natural caves, and therefore included 


in the Register of the Natural Cavities of the Italian Speleological Society. 
Anomalous cavities, constructed above the surface 
These are works built in elevation, or as part of buildings at the surface, but with characteristics similar to 


real underground spaces. 
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The best known example is documented inside the 
pyramid of Cheops at Giza (Cairo, Egypt), dating 
from the fourth dynasty, about 2500 BC [17]. The 
“Great Pyramid”, a massive structure composed of 
gigantic blocks of stone, is crossed by a series of 
tunnels and rooms where it is possible to transit for a 
total extent of about 400 metres. The oldest tunnels 
were dug under the basement, directly into the solid 
rock. They fall into the category of “underground 
cavities”. Other passages are placed in the body itself 
of the structure, above the natural level of the plateau. 
They were realized simultaneously to the laying of 
the blocks, and are therefore classified as “anomalous 
cavitie” (Fig. 12). 

Another striking example is the "Ponte Monu- 
mentale" (Monumental Bridge) at Genoa (Italy). It is 


plateau 


Fig. 12. Section of Cheops’ pyramid with tunnels and 
sepulchral chambers representation (drawing R. Bixio) 


a viaduct that crosses a road below, therefore suspended at 20 m from the ground (Fig. 14). The internal 
structure is hollow, supported by transversal wings (wall sections) in exposed stone with sub-circular arches 
that appear concentric due to the curvature of the extrados (Fig. 13). Access and investigation inside the bridge 


are faced as real speleological explorations [18]. 


3 


Fig. 13. 
Results and Discussion 


The results of the study and the indicated 
classification (Table) will enable builders, engineers, 
and urban planners to take appropriate measures to 
ensure the safety and stability of their designs already 
at the planning stage, to reduce possible risks while 
designing roads and hydro-engineering structures in 
places where underground structures are numerous. 
These investigations particularly refer to the central 
and southeastern part of Turkey, including 
Cappadocia, Akhlat, Ani, and Armenia's Syunik, 
Shirak, and Aragatsotn regions. 
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Fig. 14. Genoa (Italy). Exterior of Ponte 
Monumentale (photo A. Bixio) 
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Conclusion 


Definition of the different categories of artificial cavities according to the construction techniques, here 
illustrated, started up - as mentioned - in CNCA after decades of explorations conducted with a 
multidisciplinary approach, in various underground works all over the world, and from the need to establish 
clear and shared criteria useful to their study and classification. 

The formalized categories allow, together with the functional investigation of artificial cavities and the 
relating typological cataloguing, to have a reading-key for a basic scientific analysis, from which to develop — 
case by case - the specific research on each underground structure, aimed at better understanding the 
fundamental elements such as age, original purpose, modification, and reuse over time. 

For each artificial cavity the survey and the study of the construction techniques are fundamental, 
inextricably linked to that of the work signs left during the construction of the work itself, such as marks 
produced by pickaxes, hammers and chisels on the rock walls, from which one can also infer the digging 
direction. These observations and surveys can provide crucial elements of interpretation and understanding 
about the underground setting, borrowing and adapting to artificial cavities those systems and approaches 
improved in the last thirty years by the Archaeology of the Architecture [19-22]. They allow to read properly, 
according to stratigraphic principles, the peculiar characteristics of the “masonry evidences” shaped by 
subtraction of the raw material, following the excavation of rock masses whose modelling of empty spaces 
create the hypogean structure: from the general structures (floors, wall covering, pillars and roofs), to the 
specific elements (frames, capitals, scaffolding holes and plasters) forming the artificial cavity. 

In many cases it is possible to distinguish different “stratigraphic masonry unit” useful for relative 
dating of the structure, and often in connection with various construction techniques. The measurement of 
recurring elements in artificial cavities, such as niches and footholds, in addition to the possible presence of 
structural components also characteristic of the elevated architecture, such as brick and stone ashlars, can be 
used to recognize these elements as possible chronological indicators, 
Mensiochronology [23]. 

In summary, the study about the origin and the evolution of underground structures through the analysis of 
their way of execution and the use of the six general categories developed by the CNCA, based on the 
construction techniques above described (summarized in the Table together with the “type tree’’), is a basic 
survey instrument designed to achieve a thorough historical understanding of the hypogean architectural 
heritage. 


as developed again by the 


Table. Categories and Types of artificial cavities classified by Commissione Nazionale Cavita Artificiali of 
Societa Speleologica Italiana (R. Bixio, elaboration after Bixio and Galeazzi 2009: //document.speleo.it/) 


CONSTRUCTION TECHNIQUE 
CATEGORIES 
O EXCAVATED CAVITIES 


CAVITIES DUG IN THE SUBSOIL 
a-rock-cut structures 
b- underground structures 


CAVITIES CONSTRUCTED IN THE 
SUBSOIL 


c- tunnel excavation technique 
d- trench excavation tecnique 


MIXED ARTIFICIAL CAVITIES 
- structures dug to reach, extend 
or modify natural caves 


ANTHROPIZED CAVES 
- natural caves containing human 
interventions 


O NON-EXCAVATED CAVITIES 


CAVITIES OBTAINED BY RE-COVERING 
- underground volumes obtained by 
overlapping, embedding, burial 


ANOMALOUS CAVITIES CONSTRUCTED 


ABOVE THE SURFACE 
-strutures assimilable to undergrounds 
contained in buildings 
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Introduction 


The best solutions for heat supply have been the subject of numerous studies in the past, but the relevance 
of this research has increased due to contemporary social and political circumstances [1,2]. 

Since the RA and many other nations import fuel, their economies strongly depend on the stability of its 
price and supply, making efficient methods for heating apartment complexes a crucial concern. 

Central boiler houses and thermal power plants (TPPs), were the primary source of heat for the cities of the 
former Soviet Union and many other Eastern European nations [3]. The heat was transferred through heating 
networks. After the energy crisis of the early 1990s, when gas supplies to the country fully ceased, Armenia's 
district heating systems essentially stopped functioning. New issues with the heating of residential areas 
emerged once the gas supply was resumed in the early 2000s [4]. 

Using cutting-edge cogeneration technology at the time, some regions of the RA made an effort to 
rehabilitate central heating systems [5,6]. Due to the outdoor heating network's prolonged idleness, it 
deteriorated, and the heating systems of apartment buildings were either fully or partially disassembled. Since 
the beginning of the certificate of property rights to real estate, it has been technically impossible to restore the 
unified heat supply system for buildings because it would require large investments. It became common to 
install individual gas boilers to heat each apartment in multi-unit residential buildings as well as individual 
homes [7]. 

It is necessary to compare and assess the viability of individual and centralized heat supply options in 
apartment buildings that are currently under construction, taking into account the proportion of thermal energy 
consumed by residential buildings in the Republic of Armenia's overall energy balance as well as a noticeable 
increase in newly built residential buildings. The probability of gas equipment mishaps and water leaks in the 
heating systems is extremely high, despite the fact that individual heating systems are not supervised by the 
building operator. Due to the fact that the heating boiler, heating appliances, as well as the acquisition, 
installation, and ongoing maintenance of the system, are typically handled in this situation, the majority of 
developers continue to prefer the individual-apartment version of building heating and hot water supply [6]. 
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The theoretical value of the generated thermal energy should be utilized as the main criterion for orientation, 
however, considering the regulatory framework of the republic and the socioeconomic position of the 
population [8,9]. 

Based on the peculiarities of settlement development, individual gas boilers with a capacity of 24-32 kW 
are now used in the Republic of Armenia to heat apartments. The typical heat load of apartments or small 
private homes under "Yerevan parameters” throughout the heating season calculation period is 4-12 kW. 

According to our past research: 

e 8-12 kW in January, 
e 6-8 kW in December, 
e 8-12 kW in February, 
e 4-6 kW around in November and March. 
So, an apartment or private home's gas 
boiler, therefore, runs at 35 to 55 percent 
load [5]. The "Yerevan parameters" do not 
only apply to the city of Yerevan, it should 
be mentioned, based on the Republic of 
Armenia's territory's climate zoning, which 
takes into account climate, degree-days, the 
length of the heating season, and a seasonal 
benchmark for energy production! (Fig.1). 
The trends found by the study can be applied 
to other communities in the third zone, 
which composes around 20% of the 


republic's land area and is also the most 


densely populated in Armenia (56 percent of 1-Cold Dry, aes pees . ae 2toid hones. 
the total population) [9]. 5 - Average cold, 6 - Hot medium humid 

The choice of heat source is not necessarily ideal in a traditional district heating system with huge district 
heating boilers or TPPs because of the network's high heat loss, particularly in the case 
of TPPs. 


Since investments for medium and small businesses are mainly made through banks, attracting private 
capital is counterproductive. As a result, local residents are often forced to use individual heating systems for 
their homes. 

The main focus is on the use of individual heating and hot water systems, with the expanded use of 
renewable energy sources, in accordance with the regulation on heat supply of the program ensuring the growth 
of the energy sector of the RA until 20407”. 


Materials and Methods 


Various new building projects were surveyed in Yerevan's administrative districts of Shengavit, Kanaker- 
Zeytun, Malatia-Sebastia, Ajapnyak, and Kentron. They were contrasted in terms of how they delivered heat. 
To compare the economic viability of two systems, centralized and individual heat supply, it is crucial to 
calculate the precise values of thermal energy generated. In this study, the boiler rooms designed for multi- 
apartment buildings' own needs are considered central heating systems. 


'RACN II-7.01-2011, Construction Norms, Construction climatology, Yerevan, 2013. 

> Energy Law of the Republic of Armenia, March 7, 2001. 

3 Government of RA. Resolution No. 48-L. To adopt the Republic of Armenia's strategic plan for the development 
of the energy sector (until 2040), the plan-time-table assuring the implementation of the strategic plan for the 
development of the energy sector (until 2040), and a number of government decisions on repeal; 14.01.2021. 
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Table 1 presents the geometric and thermal characteristics of apartment buildings designed for several 
districts of Yerevan, RA in 2020-2022. The Table shows the number of stories of buildings, the number of 
apartments, the estimated number of residents, the heated volume, the calculated heating and ventilation load, 
and the specific heat and power characteristics of heating and ventilation [1]. Calculations were carried out for 
the difference in air temperatures in the room and outside (AT = 39K). The thermal protection characteristics 
of building envelopes designed after 2016 practically do not differ from each other and comply with the 
requirements of RASN 24-01-2016. When calculating the thermal and ventilation loads of buildings, the 
materials of the external structures of buildings and their design features were considered’. 


Table 1. Geometric and thermal characteristics of apartment buildings 
designed for several districts of Yerevan, RA in 2020-2022 


N 
n n S > 
- oa = < S Zy ay 
g a = 2 s 3 8 an 
@ E g E 52 3 a2X 
a oe 0 = =a Ss a af &°5 ie 
District ° 3 Ss os ne a 5 = ae 
5) on > io) io) = us} a Ss gN 
< S = o i) ° 3 Oo & = 
g 5 2 3 3 = ie eat 
5 4 5 E El = S 5° 
an 5 5 o on o 
Z Z, = Da 
1. | Shengavit 10 5300 15900 | 110 330 259.2 627 0.418 
2. | Kanaker-Zeytun 14 6048 18144 84 252 358.2 377.0 0.506 
3. | Malatia-Sebastia 18 8586 25758 72 216 447.2 471.0 0.445 
4. | Achapnyak 17 6562 19686 | 102 306 447.7 667.1 0.583 
5. | Kentron 16 11376 | 34128 | 128 384 513.4 729.6 0.386 


Let's assume that the lower combustion heat of natural gas imported to Armenia has an average value of 
8000-8200 kcal/m?. Assuming that natural gas's average lower combustion heat is 8000 kcal/m?, it turns out 
that burning 1 m? of gas results in the production of 9.3 kWh of heat. 

The studies were conducted between May 2021 and June 2022. The calculations were made taking into 
account the natural gas tariffs that were valid until 04/01/2022°, and set after 01/04/2022°. In the first case, 
with the consumption of up to 10.000 m? of natural gas, the monthly tariff for 1.000 m3 of natural gas was 
139.000 drams including value added tax, in the second case — 143. 

700 drams including value added tax. For consumers with a monthly consumption of 10.000 m3 or more, 
the price of each 1.000 m* of natural gas sold is $255.91 and $265.81, respectively, including value-added tax. 
The tariff has increased by 3-4 percent. 

During this time, the cost of natural gas sold to consumers has fluctuated significantly, as has the exchange 
rate of the US dollar against the Armenian dram (the US currency has depreciated by more than 27 percent, 
see Fig. 2, line 1). Since RA is a fuel-importing country, changes in the exchange rate have a large impact on 
the economy, especially in the gas and heat supply sectors. Taking into account exceptional fluctuations in the 
exchange rate of the US dollar (USD) against the Armenian dram (AMD) (see Fig. 2, line 2). The average 
exchange rate index for the period from May 2021 to June 202278 was also determined [10]. 


4RACN 24-01-2016, Construction Norms, Thermal Protection of Buildings, Yerevan, 2016. 

> Tigran Gnuni. 2018. Energy Balance of the Republic of Armenia. Development of Armenia’s Fourth National 
Communication and Second Biennial Update Report to the UNFCCC. 

® RACN II-7.02-1995, Construction Norms, Building Thermophysics of Fencing constructions, Yerevan, 1995. 

7 Central Bank of the Republic of Armenia, https://www.cba.am/am/SitePages/ExchangeArchive.aspx-URL 26.05.2022 

8 RA Public Services Regulatory Commission. Decision No. 83-N. 
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Fig. 2. May 2021 to June 2022 fluctuations in the USD/AMD exchange rate 


Natural gas usage for residential needs (heating, domestic hot water, cooking) is limited to 10.000 m*, while 
power plants, production facilities, boiler houses, and other facilities all use 10.000 m? or more. 

To ensure that the calculations were practical, they were also performed using the following formula 
for determining the unit cost of the product, which is currently widely used [8,9]: 


a Ki + ki + Mi + Hi 
i a, , 
K; is the annual component of capital investments (AMD/USD), 
E; is the system operation annual cost (AMD/USD), 
M, are mandatory costs regardless of production processes (such as salaries and benefits), (AMD/USD), 
H; are taxes and fees (AMD / USD), 
Q; is the annual amount of heat received (kWh/year). 


For clients with a monthly usage of up to 10.000 m3 of natural gas, the price per 1000 m? was 255.91 USD 
in 2021 and 265.85 USD in 2022”. Production facilities, massive boilers, thermal power plants, etc. use more 
than 10.000 m3 '°. 

For the observed apartment buildings, the specific value of thermal energy was calculated in 2 variants: in 
one case with a centralized heating system, in the other case with individual boilers. Table 2 shows the 
calculation results. According to studies, a medium-capacity boiler's average seasonal efficiency is about 0.85. 
The specific average cost of heat energy with centralized heat supply was around 0.042 USD/kWh (before 
01/04/2022) and around 0.044 USD/kWh (after 01/04/2022), while with decentralized heat supply, the average 
value was around 0.113 USD/kWh (until 04/01/2022) and 0.126 USD/kWh (after 04/01/2022). 


° RA Public Services Regulatory Commission. Decision N95 and, Natural Gas Supply and Use Rules, 08.07.2005. 

'0 RA Public Services Regulatory Commission. Decision Ne 221-N. - On repealing the decision 333 of November 
25, 2016 of the Public Services Regulatory Commission of the Republic of Armenia to set the tariffs for natural 
gas sold to consumers by "Gazprom Armenia" Closed Joint-Stock Company. 19.06.2020. 
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Table 2. Value of the product or specific heat capacity (SHC), (USD/ kWh) 


Value Centralized heating Individual 
Ne 1 Ne 2 Ne 3 Neo 4 Ne 5 heating 
Up to 01.04.2022 
Value of the product (SHC), 0.046 0.037 0.04 0.043 0.042 0.1 
(USD/ kWh) After 01.04.2022 
0.048 0.039 0.042 0.046 0.045 0.13 


The estimates take inflation into account in addition to variations in the dram’'s value relative to the US 
dollar and natural gas tariffs. In particular, the 12-month average inflation rate in Armenia (June 2022 
compared to June 2021) was 10.3 percent based on data released by the RA Statistical Committee. According 
to a thorough review of all indications, the cost of 1 kWh of thermal energy increased by approximately 5% 
with centralized heat supply and by approximately 12% with individual heat supply'!. 


Result and Discussion 


Based on the study, an attempt was made to justify the choice of an effective method of heat supply. In this 
paper, the authors did not consider all aspects of energy saving in buildings, which will help reduce the energy 
consumption of buildings and dependence on imported fuel and energy resources. However, having assessed 
the danger of the method of providing heat supply by individual boilers, which is widely used in multi- 
apartment buildings being designed and newly built today, attention was also paid to the energy-economic side 
of this important issue. The study carried out in the article can serve as a basis for further correction of the 
vision of heat supply in the Republic of Armenia, quickly overcoming the problems associated with the import 
of fuel resources, in particular natural gas, to create prerequisites for effective and sustainable harmonious 
development. 


Conclusion 


Thus, along with the growth of modern construction in the Republic of Armenia, important issues that 
contribute to increasing the energy independence of the country and protecting the environment are ignored. 
In the Republic of Armenia, the regulatory and technical base of approaches to the design of apartment 
buildings is far behind the requirements of the time. Not only normative documents on heat supply of multi- 
apartment buildings and new residential areas are subject to revision, but also the attitude of residents to receive 
more affordable and safe heat supply. 
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changes in any given flow direction. Through computer analysis, 

velocity change patterns in various sections along the inlet transition 

area have been constructed, allowing for the determination of fluid 

velocity at any point on the section and an estimation of the length of 

the transition area. These proposed solutions provide a framework for 

accurately constructing individual units of hydromechanical 

equipment. 

Keywords: viscous fluid, inlet section, velocity, velocity distribution, 

pressure, length. 


E-mail: asarukhanyan5 1 @mail.ru 


Introduction 


Studies of velocity distribution patterns after the inlet section in closed beds show that the particles near the 
static walls perform a decelerated motion, and the particles near the axis experience an acceleration and 
perform an accelerated motion. These two tendencies result in a rearrangement of the velocity field at a certain 
point of the movement, causing a change in the distribution pattern of velocities in the inlet section. Within the 
transition zone, the distribution of the inlet section velocities is rearranged to match the velocity distribution 
pattern of the stabilized zone in closed beds. The accuracy of the research for the transition zone of the inlet 
section is determined by the challenges associated with the precise construction of the fluid channels of the 
machinery. It is necessary to ensure clear and stable operation of the control and regulation systems. 

Therefore, the problem discussed is both relevant and of significant practical importance. The vital problem 
of fluid flow research is the mathematical model construction of the given physical phenomenon, which 
determines the applicability limits of the selected calculation method. It is vital that the built model more 
accurately describe the ongoing hydromechanical phenomena and, meanwhile provide the possibility of 
obtaining analytical solutions. 


Literature review and problem statement 


Research on the flow patterns of viscous fluid in transition areas of liquid channels, such as inlet and outlet 
sections, sudden expansion and narrowing, etc., is an applied problem that provides solutions for designing 
various machine tools. Due to their urgent and applied importance, the investigation of these problems remains 
highly relevant. 
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The starting equation for the study is the system of Navier-Stokes equations in a non-deformable medium. 
For each problem, boundary conditions are defined and the basic equations are simplified. However, obtaining 
analytical solutions to the problem is often not feasible. In such cases, modern computing techniques can 
provide an opportunity to obtain approximate solutions, which ensure accuracy in practice. 

Many theoretical and approximate calculation methods have been developed for studying the hydrodynamic 
phenomena in the transition area of the inlet section. 

Each calculation method is based on conclusions about the nature of the flow, which are used to conduct 
theoretical research and summarize the results. These conclusions often refer to certain ranges of motion, due 
to which the applicability of the obtained results is limited. A method for calculating the entrance transition 
area of plane-parallel motion was proposed in work [2], in which the author numerically implemented the 
integration of the boundary layer equations, which makes it possible to calculate the image of the distribution 
of velocities in any section of the transition area of the entrance section. Similar solutions were also 
implemented in [1]. In that study, solutions were obtained for velocity distribution and pressure change in the 
transition zone of the inlet section of the plane-parallel motion. These solutions were obtained using the 
parabolic law of velocity distribution in the boundary layer and the condition of their constancy at the core. 
Studies have been conducted to identify patterns of changes in the hydrodynamic parameters of the flow in the 
transition zone of the inlet section of the plane-parallel motion based on the approximation of the Navier- 
Stokes equations [3,4]. In these studies, a boundary problem was developed and analytical solutions were 
obtained to determine the patterns of changes in velocities and pressure. The results of the analytical solutions 
were compared with those obtained from experimental research, and the comparative analysis confirmed the 
reliability of the obtained results. 

Similar research was conducted in [9], where approximate solutions were obtained for a cylindrical and 
plane-parallel isothermal laminar motion, providing results with sufficient accuracy. The study also 
investigated the laws of the velocity field and pressure change under conditions of a four-step change of 
velocities in the boundary layer of the entrance region [14]. 

In [11], the viscous fluid flow in the inlet transition zone under the conditions of changing pressure gradient 
in axisymmetric pipes is considered. The proposed solutions are also applicable to Newtonian fluids. 

Under the conditions of varying the kinematic coefficient of viscosity according to an arbitrary law, the 
proposed solutions to the momentum equation [12] are represented by the Fourier series and the Bessel and 
Kelvin functions. Reference [13] analyzed flow stability based on the flow rate ratio to the steady-state laminar 
flow rate. This analysis led to the establishment of stability conditions for unsteady flow in a pipe. In [8], the 
authors developed a mathematical model to determine the velocity field and pressure distribution in the laminar 
motion of a viscous incompressible fluid in a two-dimensional change environment. However, the proposed 
solutions do not apply to determining the hydrodynamic parameters of the flow in areas of sudden expansion. 

The studies mentioned mainly focus on the interpretation of phenomena occurring at the entrance of the 
pipe. However, hydrodynamic parameter rearrangement phenomena also occur at other transition sites in the 
pipe, for which research is scarce. For instance, at the site of sudden expansion of the cut (D/d = 4), numerical 
integration was utilized to construct current lines using the equations of plastic fluid current movement. This 
allowed for the determination of velocity and pressure changes in the axial direction [15]. 

Extensive experimental studies on the area of sudden expansion have been conducted [16]. The authors 
utilized the magnetic-resonance tomographic method to obtain quantitative estimates of the change in 
velocities in the transition area. Under conditions of sudden, symmetric and asymmetric expansion of the 
section, quantitative estimates of the Navier-Stokes equations’ members were obtained, and the resulting 
nonlinear inhomogeneous differential equations were numerically integrated [18]. The integration results were 
compared with the results of experiments. Remarkable experimental studies on the region of sudden expansion 
were conducted in [19]. To that end, test equipment was constructed, and the areas of sudden expansion were 
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tested for the cases d/D = 0.22; 0.5 and 0.85. The investigations were carried out under both Newtonian and 
non-Newtonian fluid conditions. 

The study [7] examined the patterns of changes in the hydrodynamic parameters of laminar motion in a 
viscous fluid at the transition zone of the entrance section of a cylindrical pipe with a radius R under the 
conditions of an initially arbitrary distribution of velocities. 

Under these conditions, an axisymmetric, isothermal viscous fluid flow occurs. At the pipe inlet section, 
the velocity of the fluid u = g(r), which is moving with the velocity profile on the pipe wall, becomes zero. 


It leads to a deformation of the velocity profile, which extends over a certain distance along the length of the 
pipe. A boundary layer is formed near the pipe walls, where the velocity gradient mn becomes very large. It 
n 


causes friction forces to take on high values regardless of the viscosity coefficient i. The boundary layer 


gradually spreads from near the pipe walls until it covers the entire pipe. Therefore, studies should be 
conducted in the transition area using boundary layer equations. 


The aim and objectives of the study 


The theoretical aim of studying the transition zone is to obtain patterns of variation in hydrodynamic 
parameters along the length of the pipe and to develop methods for calculating energy losses. 

The study aims to find patterns in the hydrodynamic parameters of a viscous fluid while flowing through a 
round pipe, depending on the Reynolds number. 

To achieve this goal, the following tasks must be accomplished: 

e to develop a boundary value problem and determine the initial and boundary conditions, 

e to develop a method for solving the boundary value problem and identifying patterns of changes in 
hydrodynamic parameters of the viscous flow at the entrance section of the plane-parallel pressure 
movement, 

e toconstruct graphs of changes in axial velocities in length and Reynolds numbers, 

e to identify conditions for determining the length of the entrance section. 


Materials and Methods 
Choosing a Calculation Scheme 


The study observes the plane-parallel laminar flow of the viscous fluid at the transition zone of the entrance 
section. The origin of the Z-axis is the center of the entrance section (as shown in Fig. 1) and is directed 
infinitely long in the direction of motion. 

Plane-parallel motion in a round pipe will be considered in Cartesian coordinates, starting from the zero 
point (Fig. 1). 

The study assumes that the velocities change according to an arbitrary law at z=1 in the entrance section. It 
is necessary to find patterns of change in the hydrodynamic parameters of a viscous fluid in the transition zone, 
considering them to be axisymmetric and non-stationary. Mass forces are neglected. 


V=9 (y) “A 


Fig. 1. On the study of a viscous incompressible fluid at the inlet plane-parallel motion 
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Statement of the problem and formulation of the system of differential equations for the study 


Let's assume we have two stationary walls placed at a distance of 2 h from each other, between which there 
is a laminar, static viscous fluid flow. The pattern of the distribution of forces in the entrance section is given 


in the form of an arbitrary function V. = mal y). 
Under these conditions, an axisymmetric, isothermal movement of a viscous fluid flow occurs. Fluid 


velocity at the pipe inlet is u = p( y), and a deformation of the velocity profile occurs over a certain distance. 


A boundary layer occurs near the tube walls, where the velocity gradient a becomes very large, due to 
n 


which the friction forces assume very high values regardless of the viscosity coefficient 4 . The boundary 


layer gradually spreads from near the stationary walls and covers the entire pipe. Therefore, it is necessary to 
conduct studies in the transition zone using boundary layer equations. 

Prandtl suggested using the Nave-Stokes equations [1] for the boundary layer, which can be simplified to 
obtain the boundary layer equations. Since the main influencing forces in the boundary layer are the viscous 
forces, Prandtl simplifies the Nave-Stokes equations by ignoring terms that are very small compared to the 
viscous forces. These results in simplified equations for the boundary layer. 

Let's take the midpoint of the entrance section as the origin of the coordinates and point the oz axis in the 
direction of movement (as shown in Fig. 1). In this case, the boundary layer equations proposed by L. Prandtl 
[2] will take the following form: 


OV. av. 1 PP, dv. 


YS eV ne ere (1) 
Oz =~ sCOy p & Oy 
ov. + Ov =0. (2) 
oz (Oy 
By linearizing this system of equations according to [3,4], we will have: 
U OV...» OP ‘i Ov. (3) 


0 (= 4 ’ 
Oz po oy’ 


where U,, is the characteristic velocity of the entrance section, which is equal to the average velocity of the 


entrance section: 
U,=— | oly)dy. (5) 


Accepting that the pressures at each point of the hydraulic section have the same values, we will have: 
P=P(z). 


To integrate equations (3) and (4), the boundary conditions of the problem are defined: 


V.=0, V.=0, when y=th, (6) 
Vi= Ay), when z=0, —h<y<ch, (7) 
VV’, when z>0, —h<y<+h, (8) 


where V' is the velocity at the stabilized area, which is determined by the following equation: 
1 OP: . ‘. o’V' 
pe dy? 


(9) 
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The study assumes that the boundary layers forming near the stationary walls in the stabilized section merge 
to form a symmetric parabolic velocity profile. The solution to equation (9) will be: 


3 2 
V'(y)= u {1-2} (10) 


This corresponds to the uniform distribution of velocities in the case of plane-parallel motion [1]. 
To obtain appropriate solutions for the changes in hydrodynamic parameters in the transition zone of the 
entrance section, let's introduce dimensionless variables: 
— Vz ‘6 Zz 


Vi= 7x ,O= , Plo 
UU, h h Po 


Under these conditions, equation (3) will take the following form: 


aV(x.0)_ Po lo) 1 2°V. (12) 
aa pu; da Re a 


where Re = —2~ - is the Reynolds number. 
V 


The boundary conditions for integrating equations (12) will be: 


V..)=0,7,(«0)=y(e) EF 50, when e920, =O. 7 (0) 970). a9 
oO 


Results of research to identify patterns of changes in hydrodynamic parameters 


The solution of the inhomogeneous equation (12) under the boundary conditions (13), let's look for a 
solution in the form of a sum [5]: 


V(x,c), =U(x,0)+y(o), (14) 
where U (x, oc) represents the general solution of the homogeneous equation for inhomogeneous boundary 


conditions, while w (c) is the particular solution of the inhomogeneous equation (12) for zero boundary 


conditions. 


OU (x, 1 O°U(x, 
Let's look for the general solution of the homogeneous equation ic o) = R ie o) in the form 
Oo e€ X 


of a sum: 
U(x,o)= SC, (c)cos(y, x). (15) 


If the function U (x, oc) satisfies the equation (15), then the coefficients C, (c) must satisfy the following 


equation: 
1 
C,(o)=-—-C, (0), (16) 


from where 
y 2 
c.(0)=¢, en Zo (17) 


where C’ is the arbitrary constant coefficient. 
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By inserting the value of function C . (c) into equation (15) and taking into account equation (14), we 
obtain the following expression: 
20 2 
V(x, o)= ye C: conyarbonl Zia] swio) (18) 
k=l € 
To determine the value of the unknown function yw (c), we multiply both sides of equation (4) by bdx 


and integrate x over the interval —1< x <+1. Taking into account that the viscous liquid on the static wall 


has zero velocity, we obtain the following expression: 


ee ges (19) 


By substituting the value of velocity V, (ea co) from equation (18) into equation (19) for determining the 


function w(c), we obtain the following expression: 


2< . % 
2y'(o)= — C, siny, -exp| - a |. 20 
y'(o) Re 78s SID, i ~ (20) 
Integrating the last equation according to O , we will have: 
Ms oe Vi 
w(c)= -)'— siny, -exp| ——o |+ C(x), (21) 
kal VE Re 


where C (x) is the integration constant. 


By substituting the value of function y (c) into equation (14), we obtain the following equation for the 


velocity V, (x, a): 


00 . 2 
V.(x,0)=>°C, ody.)- sma } ox mS | +C(x). (22) 
= Y € 


k k 


= 


We determine the value of the constant C (x) from the boundary condition of the problem when 0 —> 00 
= 3(, y° _. 
Vz. (x, co) > 5 1- . The value of the constant C (x) will correspond to this boundary condition: 


we 


C(x)= =( aig?) (23) 


Thus, the final solution to the problem will be: 


roe) s 2 
V.(x,0)= ( — x?) + Sc,] cons) a fex{- Ka) (24) 
k=l 


k 


This velocity distribution equation must satisfy the boundary condition of the problem (13), according to 


which V. (x,0) = y(c). Taking into account the equation (24) and the boundary condition, we will have: 


v(o)==(-)+ Oc, conf) (25) 
k=1 


k 


To determine the values of the C,, coefficients, we must multiply both parts of the equation (25) by given 


ny, 


Si 
expression coy x)- and integrate it over the range of (-1,+1) and we will have: 


a 
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1 : 1 : 
Jaf cov, x)= srs i — > f( —x° {cow x) srs i = 
y 4 


n 


; : : (26) 
cour. — srs | coy, x)- ate la 


Vk 


; ; 0, 
J (coups) Az coups) 0H Josef ee oo Q7) 


we will have: 


Ms + 
A 
ere 


n 


As 


Lo 32) 
C, =—+—_ - —_. (28) 
sin°y, sim’ y, 


where 
1 
1) = [w(xXcos(y,x)—cosly,, )Jdx, 29) 
-l 
1 
L?) = | (1 —x* \cos(y,x)- cos(y, ))\dx , (30) 
0 


siny, 


VY; is the cos(y WJ = 0), or the positive roots of the equation 1287, = Y; . Taking into account 


k 
(28-30), we will have the solution to the problem. 


52 3 2 | pl) 372) : . 
Vsa)=3(t-s)+5| Be cosy, x) “EXP =e meee) 
= k k 


The resulting equation will satisfy the boundary conditions of the problem. 
We get the pattern of pressure change from equations (12) and (31), and we will have: 


2 j= Py 


00 Dy Re siny, 


By integrating the last equation, we can obtain the pattern of pressure variation in the transition area of the 


inlet section: 

2 © Lo (2) 2 2 

Blo) = p(0)- 22 pnt ~3t) | cxf if “| a (33) 
Po tA OY SINY, 


Let’s consider two special cases: 


1. The incoming fluid has a uniformly distributed velocity field V. (y, z)= Ay)= U,. Under these 


conditions, we will have V. (x,0) = y(x) =1, corresponding to which we will get: 


= Julayeod,x)- cosy, ))\dx = 2| (cos(y,,x)—cos(y,, ))dx = 0, 


(i- x \cos( (y,,.x)—cos( MWigea wie. 
37; 


12) = 


seat: 
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Inserting these values into the equation (31): 


V.(x0)=S(-)+ >) 


2. 
= |cos(y,x)-cos(y,))-exp] —24o |. 34) 
kal| ¥, SY, Re 


In this case, the irregularity of pressure distribution according to (33) will be equal to: 
7 _ Use 2 : ap; 
p(c) = p(0) ely =f 1 cxf Vk g Prva, (35) 
Po KAY k 
2. Let’s imagine the velocity of the liquid entering the reservoir changes according to the parabolic 


law V. (y, z) = Ay) = Al =) ), according to which we will have V, (x,0) = y(x) = (1 —x? ), 


Zz 


: . 
po) = [=x cos(y,,.x)—cos(y, Vd ican (36) 
oe 3Y; 
12) = ay (37) 
37; 


In this case, the velocity change patterns according to (31), (36) and (37) will be: 


V(x,0)= =(- x) 4 SABA) (cosy, 2)-cosy, Yen Zia. (38) 


kel OY, SINY;, 
Having the velocity distribution patterns, we can get the pressure distribution pattern in the transition zone: 
= ” U; G23-A : 3pU , 
p(c) = p(0) sel SY ( ; ), l—exp Veg da vom (39) 
Po ta 3% Re Py Re 


The resulting patterns of velocity and pressure changes make it possible to fully reveal the physical nature 
of the ongoing phenomena and make generalizations. 


Results and Discussion 


The integration of differential equations of viscous fluid flow yielded patterns of changes in the distribution 


of axial velocities V. (x, o). To visually represent these patterns of changes in the axial velocity V(x, oc) 


Zz 


across the cross-section and along the length of the transition zone, depending on the initial distribution of 


velocities V_ (x,0) = w(x) and Reynolds number, their graphs of change were constructed. Figs. 2...5 show 


Zz 


the specified calculated graphs for cases V. (x,0) =1 and V. (x,0) = (1 —x? ). 


Zz 


tite hey 


“10.01 0.02 0.04 0.06 0.08 0.2 - 
e 
Fig. 2. Graphs of changes in axial velocities V. (x, o)along the cross-section in the transition 
zone of the entrance section of the plane-parallel pressure flow at V. (x,0) =1 and 
1.2 =0.01; 2.2 =0.02; 32 = 0.04, 4.2 =0.06, 5.2 =0.08; 6 2 =0.2 
e Re Re Re Re Re 
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u/ug 0 
1.4 ‘0.1 
\0.3 
12 0.5 


=|< 


0.125 0.25 0.375 05 & 
Re 


Fig. 3. Graphs of changes in axial velocities V. (x, co) at V. (x,0) =1 and 
Lx«=0; 2.2%=0.1; 3.x=0.3; 4.x=05; 5.x=0.7; 6.x=0.9 


-1 0.08 0.1 


Fig. 4. Graphs of changes of the axial velocities V. (x, co) along the cross-section on the 


entrance section of the plane-parallel pressure flow at V. (x,0) = (1 = x’) and 


1.2 =0.08; 2.2 =0.1; 3.2 =0.2; 4.2 =05 
Pe Pe Pe Pe 


5 & 
0.05 Re 


Fig. 5. Graphs of changes in axial velocities V. (x, o) at V. (x,0) = (1 — x?) and 
1.x=0; 2.x=0.1; 3.x=0.3; 4.x=0.5; 5.x=0.7; 6.x=0.9. 
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The dynamics of the transition zone were determined based on the results of numerical calculations and 
graphs, using the condition that the ratio of axial velocities between the transitional and stabilized areas is 0.99. 


From this condition, it was found that L = =e =0.147. The proposed method for calculating velocity 
e€ 


rearrangements in the transition area of the entrance section enables the determination of patterns of changes 
in hydrodynamic parameters of the flow under general boundary conditions. Using the identified relationships, 
the deformation process of the velocity field at the transition site of the entrance section was determined for 
both constant and parabolic law distributions of incoming liquid. This enables the calculation of changes in 
the hydrodynamic parameters of the flow and facilitates generalizations of the results. 

The patterns of changes in the velocity field in the transition area of the entrance section of plane-parallel 
pressure flow, along with the graphs constructed based on these findings, provide valuable information for the 
accurate design of hydromechanical equipment units. 


Conclusion 


The study of viscous fluid flow in the transition zone of the entrance area of a plane-parallel pressure flow 
has revealed the following findings: 

e A boundary value problem was formulated to investigate the regularities of changes in hydrodynamic 
parameters of a viscous incompressible fluid in the transition zone of the entrance section of a plane- 
parallel pressure flow. 

e A method for solving the boundary value problem was developed, and formulas were derived to calculate 
axial velocities across the cross-section and pressure along the length of the diffuser. 

e Graphs were constructed to illustrate the changes in hydrodynamic parameters of the flow across the cross- 
section and along the length of the transition zone for constant and parabolic distributions of initial 
velocities at the entrance section. 

e A calculation formula was derived to determine the length of the transition zone of the inlet section of a 
plane-parallel pressure flow. 

The solutions derived from the approximating Navier-Stokes equations to identify patterns of changes in 
hydrodynamic parameters of a plane-parallel pressure flow with constant and parabolic distributions of initial 
velocities at the entrance section provide a comprehensive understanding of the underlying processes. These 
results are essential for hydraulic calculations of various systems. 
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Abstract: This article presents the use of modified pumice as an Marine Kalantaryan* 
environmentally friendly adsorbent for copper (II) removal from E-mail: kalantaryanm @ mail.ru 
wastewater. The water pollution by toxic elements is a major concern 

for human health and environmental quality. New and cheaper 


methods of wastewater treatment are increasing the quality of the Received: 02.04.2023 
environment and reducing negative impacts on fauna, flora, and Revised: 20.04.2023 
human beings. The sorption technique is considered a cost-effective Accepted: 02.05.2023 


method for effectively removing heavy metals. In recent years, there 
have been increasing studies dedicated to using low-cost adsorbents 
such as pumice. For the study, Kuchak pumice has been used. The © The Author(s) 2023 
modified pumice was prepared by surface modification with 
polysiloxane, evaluated by studying the effects of pH, contact time, 
dosage, and initial concentration, and was optimized in batch 


processing mode. The chemical changes in pumice were fully This work is licensed under a Creative 
characterized using FT-IR techniques. Overall, these results suggest Commons Attribution-NonCommercial 
that surface-modified pumice is a low-cost adsorbent for the removal 4.0 International License 

of copper (II). 


Keywords: pumice, modification, polysiloxane, adsorbent, heavy 
metal removal. 


Introduction 


The rapid development of the industry has resulted in a significant increase in heavy metal contamination, 
which has become a major environmental issue. The most significant challenge facing humanity today is 
effectively managing the problem of water pollution caused by heavy metals. The physical and chemical 
techniques traditionally used for cleanup are associated with high capital costs and can damage already 
contaminated areas. In today's highly industrialized society, several industrial operations release aqueous 
effluents containing heavy metals. The presence of heavy metals in aquatic systems can harm numerous living 
species [1]. Copper (Cu**) ions are present in many biological systems at low levels but can be toxic at high 
concentrations. Although copper can exist in different forms, such as Cu (0), Cu (1), and Cu (II), Cu (ID) is the 
primary species of concern in aqueous solutions. Cu* rapidly binds to organic and inorganic compounds in 
aqueous solutions, which is pH dependent. It harms aquatic life and makes natural water unsuitable for human 
use. Ingesting high concentrations of Cu can also be toxic to humans, leading to cancer and promoting 
oxidation [2,3]. Toxic metals can now be removed from soil and water using various clean-up technologies 
developed over time. Currently, physical-chemical processes such as filtration, chemical precipitation, ion 
exchange, adsorption, and electrodeposition, form the foundation of the most extensively utilized remediation 
technologies [3]. Different types of materials have been used as adsorbents for copper adsorption, such as 
activated carbons and zeolites [4], alumina - silica [5], molecular sieve powder [6] and different pumices. The 
current study aims to establish the ideal conditions for the maximum amount of copper to be absorbed by 
pumice from aqueous solution and to evaluate the effects of various experimental parameters on copper 
adsorption. 

Different treatment strategies, such as physical, chemical, and biological treatments, were discussed based 
on the past few years. These introduced methods are adsorption, membrane filtration, electrodialysis, and 
photocatalysis. However, there are some precautions that should be taken seriously as they are influenced by 
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parameters such as the initial concentration of copper ions, pH values, economic parameters such as operation 
cost, and the environmental effects and compatibility of each of the various methods conducted. 

Currently, many adsorbents are being produced and used for cleaning up wastewater. These adsorbents can 
be of either natural or artificial origin. Considering the numerous advantages of natural materials, including 
their availability in large quantities and non-toxicity, it is more appropriate to use adsorbents based on natural 
materials. Efforts have been made recently to modify the surface of these natural adsorbents to increase their 
efficiency, and extensive research has been conducted in this direction [1,2]. 

The aim of the present study is to investigate the influence of various experimental parameters on copper 
adsorption and determine the optimum conditions for maximum adsorption of copper by pumice from aqueous 
solution. The influence of temperature, sorbent mass, solution pH and sorbent chemical modification on 
adsorption process is discussed. 

The Republic of Armenia is rich in non-metallic mineral diversity and abundance. The light rocks (tuffs, 
perlite, pumice stone, zeolite, scoria, etc.) generated because of volcanic activities in the mountains of Armenia 
are of particular value and significance. 

In this article Kuchak pumice deposit is used as an adsorbent for copper removal from aqueous media. The 
pumice on the territory of the Republic of Armenia according to their petrographic, physical and mechanical 
characteristics are divided into two types: Ani and lithoidal Kuchak pumice is one of the Ani type pumice 
varieties. It is located at an altitude of 2050 meters above sea level. On the southern side of the pumice deposit 
there is a powerful layer of pumice grains with a capacity of 6 ...7 meters and pumicite with a capacity of 
4 ...6 meters. Studies have shown that the pumice is composed of alumicillicates in which the alkali oxide 
content varies from 1.5 to 5%, SiOz from 71% to 75%, and Al2O3 from 12% to 14%. In the Kuchak pumice 
deposit samples, the content of alkaline oxides varies from 0.03 to 0.1%, SiOz from 70% to 73% and Al203 
from 12 to 16% [7]. 


Materials and Methods 


The analysis of the data obtained shows that the Kuchak pumice deposit is environmentally safe and is 
a chemically neutral silicate rock in water medium. Kuchak pumice has a porosity of 1.64 to 32.70 um 
(Fig. 1). The total porosity is about 72.2 - 79.4 % [7]. 

The presented data show that Kuchak mine pumice is an aluminosilicate rock with well-developed porosity, 
mechanical strength, high buoyancy. It is chemically inert and eco-friendly. 

In this study the following reagents were used: polysiloxane emulsion, copper sulphate (CuSO4.5H20), and 
hydrochloric acid (HCl). All the reagents used were of highly pure grade. The deionized water was used for 
all reagent solutions. 

Polysiloxane has an extensive record of success and is frequently used in medical applications. Elastomers, 
gels, lubricants, foams, and adhesives are examples of different material types. Polysiloxane is relatively 
stable and chemically inert. They have minimal moisture uptake and are hydrophobic. They offer excellent 
electrical insulation qualities. Nearly all polysiloxane is based on polymethylsiloxane (Fig. 2). 


CH, 
| 
[-O-Si-], 
CH, 


DET: SE Detector 


50 um 
Device: TS5130mm 


Fig. 1. Pore sizes of Kuchak pumice deposit [8] Fig. 2. Chemical structure of polysiloxanes 
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Polysiloxanes (PDMS) or silicones, are a general category of polymers consisting of a silicon-oxygen 
backbone with organic groups, typically methyl groups, attached to the silicon atoms [9,10,11]. Organic side 
groups can be used to link two or more chains together. By varying the -Si-O- chain length, side groups, and 
crosslinking extent, silicones with properties ranging from liquids to hard plastics can be synthesized. Silicone 
synthesis typically involves the hydrolysis of chlorosilanes into linear or cyclic siloxane oligomers, which are 
then polymerized into polysiloxanes by polycondensation or polymerization, respectively. Polysiloxanes, 
characterized by unique material properties combining biocompatibility and biodurability, have found 
widespread application in healthcare [12,13,14,15]. 

About 1.0 g of polysiloxane was mixed with 10 g of pumice and 30 mL of distilled water. First, the pumice 
was washed with water and dried in an oven at 105°C for one hour. Then the polysiloxane emulsion was mixed 
in 1:30 ratio with distilled water to form a stable mixture. The mixture was then stirred for 5 hours at room 
temperature. Then it was filtered from the solution, washed, and dried in an oven at 70 °C for 3 hours to allow 
the emulsion to adhere to the surface. Finally, the modified pumice was washed with water to remove excess 
emulsion and dried again. Different particle sizes (2.36, 3.54, and 4.87 mm) were collected and used as an 
adsorbent for the analysis. All pH values were measured with a pH meter (HACH LANGE HQ 14d). Cu”* 
solutions of different concentrations (5 mg/L, 10 mg/L, 50 mg/L, and 100 mg/L) were prepared in deionized 
water using CuSO.. Batch adsorption equilibrium studies of Cu were carried out at room temperature in 
beakers filled with 2.5 g of modified pumice. The contact times were one or two hours. The pH dependence 
experiment was carried out as a preliminary experiment to establish if the material possesses any affinity for 
Cu (ID ions. The pH of the solution is a known factor with significant impact on metal ion removal processes. 
The experiment was conducted in the pH range 2-6. 

Fig. 3 shows the effect of pH on copper adsorption onto pumice. The adsorption capacity of pumice samples 
increased with increasing pH. The highest adsorption value is observed at a pH of four. 
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Fig. 3. Effect of pH on the adsorption of copper by surface-modified pumice 


Results and Discussion 


The effect of temperature on the adsorption of copper was studied by conducting a series of experiments at 
25, 30, and 40 °C. There is a small effect of temperature on the adsorption of copper by pumice in aqueous 
solution. Because elevated temperature would be costly in commercial applications, experiments were run at 
25 °C. The functional groups were found from the pumice and polysiloxane-modified samples using FT-IR 
spectrometer. The FT-IR spectrum of pumice is shown in Fig. 4. Compared to that of pumice, Si-CHs3 stretching 
bonds are observed for modified pumice under 1260 cm! and (Si-O-) bonds under 1100 cm’. 
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Fig. 4. FT-IR spectrum of modified and unmodified pumices (PerklinElmer, Spectrum Two) 


The siloxane groups on the modified pumice surface can form coordinate bonds with metal ions such as 
copper (Cu?*), allowing for effective removal from solution. The adsorption performance of polysiloxane- 
modified pumice will depend on various factors, including the surface area, pore size, concentration of metal 
ions, pH, and contact time. 

Copper concentrations before and after adsorption were measured using UV-Vis Spectrophotometer (Cary- 
60). The wavelength of 10.0 nm''corresponds to the maximum absorbance of the copper on pumice. 
The spectra show after modification by polysiloxane the copper absorbance increased by 65% compared with 


unmodified pumice (Fig. 5). 
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Fig. 5. UV-Vis Spectrophotometric data (Cary-60) 
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Conclusion 


The efficiency of the pumice for the adsorption of copper from an aqueous solution was investigated using 
batch adsorption technique under different experimental conditions. Reported results showed that the 
adsorption varied strongly by pH. Temperature had a minor effect on the adsorption of copper by pumice. Due 
to its low cost and ready availability, it can be used as an efficient adsorbent material for the adsorption of 
copper from contaminated aqueous solutions. The modification of pumice with PDMS introduces siloxane 
(-Si-O-) functional groups onto the surface of the pumice particles. These siloxane groups can provide 
adsorption sites for metal ions due to their ability to coordinate with metal cations. 

The high surface area and porous nature of pumice, combined with the adsorption properties of PDMS, can 
enhance the adsorption capacity for metal ions. The siloxane groups on the modified pumice surface can form 
coordinate bonds with metal ions such as copper (Cu**), allowing for effective removal from solution. The 
adsorption performance of polysiloxane modified pumice will depend on various factors, including the surface 
area, pore size, concentration of metal ions, pH, and contact time. Additionally, the specific synthesis and 
modification methods of the polysiloxane-modified pumice may influence its adsorption properties. 

Overall, polysiloxane-modified pumice has the potential to be an effective adsorbent for copper ions. 
However, the actual adsorption capacity and efficiency would require experimental evaluation and 
optimization for specific conditions and metal ion concentrations. In our future research, we aim to develop 


methods to extract copper ions from polysiloxane-modified adsorbents. 
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